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ABSTRACT 


This  document  describes  the  complete  design  of  the  Naval 
Atmospheric  Boundary  Layer  (NABL)  Model  System  in  terms  of  its 
Computer  Software  Components  (CSCs)  and  Computer  Software  Units 
(CSUs). 
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NAVY  ATMOSPHERIC  BOUNDARY  LAYER  (NABL)  MODEL  SYSTEM: 
SOFTWARE/ INTERFACE  DESIGN  DOCUMENT 


1  SCOPE. 

1.1  Identification.  This  design  document  applies  to  the  Naval 
Atmospheric  Boundary  Layer  (NABL)  Model  System  which  is  identi¬ 
fied  and  abbreviated  as  NABLS.  Its  identification  number  is  "to 
be  determined".  This  design  is  derived  from  the  NABLS  require¬ 
ments  document  (reference  a,  in  section  2.  below). 

1.2  System  overview.  The  NABLS  will  implement  the  "NABL  model" 
on  the  Fleet  Numerical  Oceanography  Center  (FNOC),  Monterey,  CA 
computer  system.  The  NABL  model  was  developed  by  the  Naval 
Oceanographic  and  Atmospheric  Research  Laboratory  (NOARL),  Atmos¬ 
pheric  Directorate,  Monterey  CA  to  produce  short-range  (0  to  36 
hour)  forecasts  of  meteorological  conditions  in  the  atmospheric 
boundary  layer.  Although  NABL  is  a  one-dimensional  "stick" 
model,  it  uses  turbulence  paramet r i zat i ons  and  a  vertical  resolu¬ 
tion  which  are  superior  to  those  in  the  current  Navy  operational 
regional  and  global  atmospheric  prediction  model  systems  (NORAPS 
and  NOGAPS  respectively).  NABL  obtains  adiabatic  tendencies  of 
temperature  and  moisture  from  NORAPS  or  NOGAPS,  giving  it  three- 
dimensional  character,  although  NABL  does  not  interact  in  a  two- 
way  sense  with  either  of  the  larger-scale  prediction  model  sys¬ 
tems.  NABL  has  demonstrated  superior  skill  at  predicting  the 
vertical  structure  of  temperature  and  refractivity  in  the  bound¬ 
ary  layer  (refs,  g  and  h.). 

The  NABLS  will  permit  a  system  user  to  specify  where  geo¬ 
graphically  and  for  what  forecast  intervals  and  times  the  NABL 
model  will  be  run;  will  cause  the  requisite  initial  condition  and 
tendency  information  from  one  of  the  larger-scale  models  to  be 
made  available  for  input  to  NABLS;  and  will  provide  for  saving 
the  resultant,  high-resolution  NABL  model  output  in  a  convenient 
format(s)  for  use  by  such  other  applications  (external  to  NABLS) 
which  may  now  or  in  the  future  have  need  for  such  information. 

1.3  Document  overview.  This  document  describes  the  complete 
design  of  the  NABLS  in  terms  of  its  Computer  Software  Components 
(CSCs)  and  Computer  Software  Units  (CSUs).  It  would  normally  be 
used  as  a  basis  for  coding  each  CSU  but  that  is  not  true  in  this 
instance  for  reasons  explained  in  the  Note  at  the  end  of  this 
subparagraph.  This  document  is  organized  according  to  the  DOD 
2167A  standard  for  the  Software  Design  Document  (SDD) 
(ref.  b,c.).  That  SRS  standard  has  been  modified  slightly  to 
incorporate  into  that  one  document  the  pertinent,  non-r edundant 
sections  of  the  DOD  2167A  standard  for  the  Interface  Design 
Document  (IDD)  (ret.  d.).  This  modification  caused  section  3  of 
the  SDD  outline  to  be  expanded  to  include  an  additional  subsec¬ 
tion  3.3  entitled  Interface  Description. 

Section  1  identifies  and  provides  an  overview  of  the  NABLS. 
Section  2  lists  applicable  documents.  Section  3  provides  an 
overview  of  the  NABLS,  describes  its  various  components,  and 
provides  a  preliminary  design  in  terms  of  execution  control  and 
data  flow  between  those  components.  Section  4  provided  design 
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details  to  the  unit  level.  Section  5  describes  the  global  data 
elements  within  the  NABLS.  Section  6  describes  the  data  files 
which  are  shared  between  NABLS  components.  Section  7  provides 
traceability  of  the  requirements  stated  in  the  Requirements 
Specification  (ref.  a)  to  the  paragraphs  in  this  document. 

Note:  At  the  time  this  document  was  being  prepared,  the 

NABL  model  and  various  interfaces  thereto  had  already  been  in¬ 
stalled  at  FNOC  and  the  "system"  as  such  was  undergoing  informal 
operational  checkout.  This  design  therefore  documents  the  re¬ 
quirements  for  and  architecture  of  the  NABLS  retrospectively. 
Accordingly,  emphasis  herein  is  given  to  providing  information 
which  will  benefit  the  maintenance  programmer  for  the  system  as 
it  now  exists  and/or  the  analyst  who  may  wish  to  make  improve¬ 
ments  to  NABLS  as  it  is  now  implemented.  Requirements  and  design 
specifications  in  this  document  which  vary  substantially  from 
NABLS  as  it  is  now  implemented  are  noted  (in  bold  face)  as  poten¬ 
tial  changes  or  enhancements  which  could  be  made  at  a  later  date. 

2.  REFERENCED  DOCUMENTS. 

2.1  Government  documents.  The  following  documents  of  the  exact 
issue  shown  form  a  part  of  this  specification  to  the  extent 
specified  herein.  In  the  event  of  conflict  between  the  documents 
referenced  herein  and  the  contents  of  this  specification,  the 
contents  of  this  specification  shall  be  considered  a  superseding 
requirement . 

a .  Navy  Atmospheric  Boundary  Laver  (NABL)  Model  System: 
Software/Interface  Requirements  Specification  (NABLS  SRS/IRS). 
NOARL  Technical  Note  182,  November  1991. 

b.  DOD-STD-2 1 67A,  Military  Standard,  Defense  System  Software 
Development .  DoD,  Washington  D.C.,  29  February  1988. 

c.  DI-MCCR-80012A,  Data  Item  Description  for  SDD. 

d.  DI -MCCR-80027A ,  Data  Item  Description  for  IDD. 

e.  FNOC  Computer  User  Guide  Edition  2  (updated). 

f.  FNOC  Cyber  205  User  Guide,  December  1989. 

g.  Vogel,  G.  N.,  and  R.  H.  Langland,  1991:  An  Evaluation  of 
the  Navy  Atmospheric  Boundary  Laver  (NABL)  Model.  NOARL  Techni¬ 
cal  Note  160,  August  1991. 

h.  Vogel,  G.  N.,  1991;  Assessment  and  Short-range  Forecast¬ 
ing  of  Tropospheric  Ducting;  Comparison  of  Results  from  the 
Navy's  Local.  Regional  and  Global  Forecast  Model.  NOARL  Report 
23,  1991. 

i.  Lewit,  H.  (CSC),  NGPOINT.  The  NABL  NOGAPS  Geographic 
Point  File  Interface  System.  User’s  Instructions.  18  Dec  1990. 
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j.  Lewit,  H.  (CSC),  XRPNABL.  The  NABL  NORAPS  Geographic 
Point  File  Interface  System.  User's  Instructions.  20  Dec  1990. 

k.  FLENUMOCEAXCENINST  5234.5  (CH-1),  FLEXUMOCEAXCEN  Software 
Standards .  12  January  1990. 

Copies  of  specifications,  standards,  drawings,  and  publications 
required  by  suppliers  in  connection  with  specified  procurement 
functions  should  be  obtained  from  the  contracting  agency  or  as 
directed  by  the  contracting  officer. 

2.2  Non-Government  documents.  The  following  documents  of  the 
exact  issue  shown  form  a  part  of  this  specification  to  the  extent 
specified  herein.  In  the  event  of  conflict  between  the  documents 
referenced  herein  and  the  contents  of  this  specification,  the 
contents  of  this  specification  shall  be  considered  a  superseding 
r equi rement . 

l.  CDC  XOS/BE  Reference  Manual  (available  from  FNOC/NOARL 
1 ibrary ) . 

m.  CDC  VSOS  Version  2  Reference  Manual,  Volume  1  (available 
from  FXOC/XOARL  library). 

n.  CDC  FORTRAN  200  Version  1  Reference  Manual  (available 
from  FXOC/XOARL  library). 
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3.  PRELIMINARY  DESIGN.  Note:  This  section  would  normally  be 
completed  and  then  reviewed  as  a  true  "preliminary"  cut  at  system 
design.  The  "detailed"  design  sections  which  follow  would  be 
based  on  that  design  review.  Because  this  document  is  being 
prepared  retrospectively,  this  section  provides  a  top-level  view 
of  the  existing  system  and  is  in  no  way  preliminary.  The  section 
title  is  retained,  however,  to  comply  with  the  documentation 
standard  ( ref .  c) . 

3.1  NABLS  overview.  Figures  3-1,  3-2  and  3-3  are  reproduced 
herein  from  ref.  a.  Figure  3-1  depicts  the  relationship  of  the 
NABLS  to  its  several  external  interfaces.  The  figure  is  also 
known  as  the  context  diagram  because  it  shows  the  domain  of  the 
system.  Detailed  specifications  for  each  external  interface  are 
provided  in  ref.  a. 

Figures  3-2  and  3-3  present  data  flow  within  the  NABLS.  The 
symbology  used  in  these  data  flow  diagrams  (DFDs)  is  described  in 
Appendix  A. 

3.1.1  NABLS  architecture.  Figure  3-4  shows  the  overall  NABLS 
architecture.  It  relates  each  requirement  (keyed  to  the  identi¬ 
fying  numbers  assigned  to  the  "bubbles"  in  figures  3-2  and  3-3) 
to  the  programs  and  procedures  which  support  those  requirements. 
Note:  r equi r ement/bubbl e  1.4  (Transfer  Model  Output)  is  not  shown 
because  of  its  loose  relationship  to  the  rest  of  NABLS,  as  dis¬ 
cussed  in  paragraph  3. 2. 1.4  of  ref.  a. 

The  programs  and  procedures  shown  on  figure  3-4  constitute 
the  top-level  Computer  Software  Components  (CSCs)  of  the  NABLS. 
They  are  highlighted  by  their  double-sided  boxes.  All  of  these 
CSCs  already  exist  on  the  FNOC  computer  system.  Some,  as  dis¬ 
cussed  in  ref.  a  and  as  will  be  noted  again  in  this  document,  are 
recommended  for  modification  and/or  consolidation. 

Note:  In  this  document  the  CSCs  are  discussed  as  separate 
entities  as  they  existed  on  or  about  25  October  1991.  That 
caveat  is  required  since,  even  as  this  document  was  being  draft¬ 
ed,  changes  were  being  made  to  program  NABL  to  make  that  primary 
element  of  NABLS  run  more  efficiently  on  both  the  FED  ("front- 
end")  and  SAM  ("super")  computer  systems.  Some  of  those  recent 
changes  were  in  direct  response  to  recommendations  noted  in 
r ef .  a . 


3.1.2  System  states  and  modes.  The  NABLS,  and  in  particular 
program  NABL,  is  designed  to  operate  in  the  batch  mode  as  part  of 
a  scheduled  job(s)  -  either  under  NOS/BE  ((CDC's)  Network 
Operating  System/Batch  Environment  (ref.  1))  on  FED  and/or  under 
VSOS  ((CDC's)  Virtual  Storage  Operating  System  (ref.  m) )  on  SAM. 
Program  NABL  satisfies  requirements/bubbles  2  through  4  in  figure 
3-2. 

The  system  may  also  be  operated  in  an  unscheduled  batch 
mode.  This  would  be  the  usual  mode  of  exercising  the  INTERFACE 
LARGER-SCALE  MODELS  capability  ( requi rement/bubbl e  1)  which  re¬ 
quires  interactive  input  from  a  user;  but,  only  when  the  user 
initiates  the  process.  Program  NABL  may  also  be  exercised  in  an 
unscheduled  batch  mode  for  R&D  applications. 
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Figure  3-3.  DFO  of  Interface  Larger-scale  Models  (Bubble  T'. 
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Figure  3-4.  The  NflBLS  ftrchi lecture. 
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The  NAMELIST_DATA  which  accompanies  the  NABL  initiating  job 
will  also  be  updated  on  an  unscheduled  basis.  This  will  be 
accomplished  by  the  responsible  program  maintenance  person  pre¬ 
paring  and  submitting  an  AUTOF  update  to  the  AUTO  File  as  de¬ 
scribed  in  ref.  e  (section  7.6). 

3.1.3  Memory  and  processing  time  allocation.  The  NABLS  CSCs 
will  be  governed  by  the  programming  standards  and  resource  usage 
provisions  of  refs,  e  (section  6)  and  f. 

3.2  NABLS  design  description.  This  section  provides  a  design 
description  of  the  NABLS  CSCs.  Each  subparagraph  below  corre¬ 
sponds  to  one  of  the  seven  programs/procedures  identified  by  a 
double-sided  box  in  figure  3-4. 

3.2.1  Program  NABL  (NABL).  This  batch-mode  atmospheric  predic¬ 
tion  program  is  central  to  the  NABLS.  This  one-dimensional 
boundary  layer  model  is  initialized  with  profiles  of  wind,  tem¬ 
perature,  and  moisture  at  preselected  locations  as  provided  from 
either  NORAPS  or  NOGAPS.  The  evolution  of  these  profiles  is 
forecast  based  on  parametrized  physical  processes,  surface  forc¬ 
ing,  and  synoptic-scale  adiabatic  advective  tendencies. 

Program  NABL  satisfies  the  Process  Input  Data,  Compute 
Forecast  Variables,  and  Prepare  Output  requirements  (paragraphs 
3.2.2,  3.2.3  and  3.2.4  of  ref.  a).  A  table  of  the  COMPLETE_FORE- 
CAST_DATA  file  data  elements  may  be  found  in  Appendix  D. 

Control  of  program  execution  is  provided  by  the  NAMELIST 
data  read  in  by  NABL  from  the  NAMEL I ST_DATA  i  n t e r f ac e / f i 1 e 
(paragraph  3.1.2  of  ref.  a).  A  table  of  the  NAMELIST_DATA  file 
data  elements  may  be  found  in  Appendix  D. 

The  individual  modules,  the  Computer  Software  Units  (CSUs), 
which  comprise  the  NABL  program  are  presented  in  Section  4.1. 

3.2.2  Procedure  NGPOINT  (NGPOINT).  This  Cyber  Control  Language 
(CCL)  procedure  is  initiated  by  a  user  who  wishes  to  input  a  list 
or  a  change  to  an  existing  list  of  NABLS  geographic  point  loca¬ 
tions  for  input  to  the  NOGAPS  synoptic  scale  atmospheric  predic¬ 
tion  system.  This  list  will  control  the  location  of  profiles 
prepared  by  NOGAPS  for  input  to  program  NABL. 

This  procedure,  which  is  described  more  fully  in  ref.  i, 
starts  the  interactive  interface  program  NGNPNT.  If  warranted, 
based  upon  "sense  switch"  setting  provided  by  that  program,  this 
procedure  will  cause  the  resultant  DATAPOINT_NAMELI ST  to  be 
transferred  into  the  SAM  "super"  computer  system  for  use  the  next 
time  the  NOGAPS  is  executed. 

Procedure  NGPOINT  satisfies  the  Interface  Control  require¬ 
ment  (paragraph  3.2. 1.1  or  ref.  a)  for  NOGAPS. 

Control  of  procedure  execution  is  provided  by  the  sense 
switch  settings  mentioned  above  -  only  if  SWl  is  set  ON  is  the 
NAMELIST  information  transfer  initiated. 

3.2.3  Procedure  NRPNABL  (NRPNABL).  This  CCL  procedure,  which  is 
described  more  fully  in  ref.  j,  is  the  NORAPS  counterpart  of 
NGPOINT  above.  Note:  as  discussed  in  ref.  a,  the  two  procedures 
should  probably  be  combined  at  some  convenient  time. 
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If/when  this  procedure  is  combined  with  procedure  NGPOINT, 
sense  switch  two  could  be  set  ON  by  the  interactive  program  to 
indicate  that  a  NORAPS  instead  of  a  NOGAPS  namelist  is  involved. 
This  switch  setting  would  then  be  used  to  initiate  the  appropri¬ 
ate  file  transfer  since  different  procedures  are  used  for  that 
purpose  -  depending  on  which  of  the  two  larger-scale  model  sys¬ 
tems  is  involved. 

3.2.4  Program  NGNPNT  (NGNPNT).  This  interactive  program  is 
initiated  by  procedure  NGPOINT  above.  By  query/response  it 
obtains  the  user’s  input  with  regard  to  which  geographic  points 
should  be  added  to,  deleted  from  or  changed  on  the  current  DATA- 
POINT_DETAILS  file.  When  the  user  has  completed  this  dialogue, 
the  programs  writes  out  a  new  DATAPOINT_NAMEL I  ST  file  (assuming 
there  have  been  changes,  or  that  it  is  an  entirely  new  file)  and 
sets  the  appropriate  sense  switch(es)  to  inform  procedure  NGPOINT 
what  has  happened. 

This  program  CATALOGS  the  new/revised  DATAPOINT_DETAI LS  file 
if  one  is  created  (the  expected  case).  A  table  of  the  DATA- 
P01NT_DETAILS  and  DATAPOINT_NAMELI ST  file  data  elements  may  be 
found  in  Appendix  D. 

Program  NGNPNT  satisfies  the  Get  User  Input  requirement 
(paragraph  3. 2. 1.2  of  ref.  a)  for  NOGAPS. 

Control  of  program  execution  is  provided  entirely  by  user 
responses . 

The  individual  CSUs  (modules)  comprising  the  NGNPNT  program 
are  presented  in  Section  4.4. 

3.2.5  Program  NRPPNT  (NRPPNT).  This  interactive  program  is  the 
NORAPS  counterpart  of  NGNPNT  above. 

This  program  differs  from  its  NOGAPS  counterpart  above  by 
writing  an  "intermediary"  CCL  procedure,  PRO,  which  will  pass  on 
the  identification  of  the  particular  NORAPS  area/domain  the  user 
has  specified.  This  is  not  necessary  in  the  case  of  NOGAPS  since 
there  is  only  one  (global)  domain.  A  second  difference  is  caused 
by  the  current  (questionable)  convention  of  not  requiring  users 
to  input  their  assessment  as  to  the  "wet"  or  "dry"  condition  of 
each  geographic  point.  A  third  difference  is  caused  by  the  need 
to  append  at  the  beginning  of  each  DATAPOINT_NA.MELIST  file  a  set 
of  20  integer,  logical,  or  character  variables  which  control 
global  (programmatic)  aspects  of  NORAPS  model  execution.  These 
variables  are  unrelated  to  NABL  and  can  vary  from  NORAPS  domain 
to  NORAPS  domain. 

The  individual  CSUs  (modules)  comprising  the  NRPPNT  program 
are  presented  in  Section  4.5. 

Note:  as  discussed  in  ref.  a,  the  two  programs  (NGNPNT  and 
NRPPNT)  should  be  combined  at  some  convenient  time.  (See  related 
Note:  in  paragraph  3. 2. 1.2  of  ref.  a  and  Appendix  E  hereto.) 

3.2.6  Procedure  NGNMOD  (NGNMOD).  This  CCL  procedure  transfers  a 
new  or  revised  list  of  geographic  points  ( DATAPOINT_NAMEL I  ST)  to 
SAM  for  use  by  NOGAPS.  It  ROUTES  a  job  (OKNMD)  into  SAM’s  input 
queue  with  the  new  data  file  attached  to  the  job.  That  job  will 
COPY  and  SWITCH  the  tag-a-long  file  to  replace  the  older  DATA- 
POINT_NAMELIST  file. 
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Procedure  NGNMOD  satisfies  the  Transfer  Points  requirement 
(paragraph  3. 2. 1.3  of  ref.  a)  for  NOGAPS. 

Control  of  procedure  NGNMOD  initiation  is  provided  by  proce¬ 
dure  NGPOINT.  Control  during  procedure  execution  is  not  re¬ 
quired. 

3.2.7  Procedure  NRNMOD  (NRNMOD).  This  CCL  procedure  is  the 
NORAPS  counterpart  of  NGNMOD  above. 

It  differs  from  NGNMOD  because  it  must  be  aware  of  which 
NORAPS  area  is  to  be  changed  and  then  pass  this  information  on 
into  SAM  so  that  the  correct  area-corresponding  namelist  file  can 
be  updated.  An  integer  area  identifier  is  passed  into  this 
"permanent”  procedure  as  a  parameter  when  this  procedure  is 
called  (BEGINed)  from  within  the  "intermediary"  CCL  procedure 
named  PRO  which  was  discussed  above  in  paragraph  3.2.5. 

Note:  these  two  permanent  transfer  procedures  (NGNMOD  and 
NRNMOD)  are  substantially  different.  There  would  be  no  benefit 
if  they  were  combined. 

3.3  NABLS  Interface  descriptions.  Figure  3-1  shows  the  external 
interfaces  of  the  NABLS.  Descriptions  of  each  of  those  inter¬ 
faces  is  provided  in  the  following  subparagraphs.  (The  informa¬ 
tion  in  this  subsection  draws  heavily  on  material  presented  in 
subsection  3.1  of  ref.  a.) 

3.3.1  USER_ONLINE_TERMINAL .  This  interface  provides  interactive 
input  from  the  user  who  is  prompted  to  provide  details  concerning 
the  location  and  character  of  the  data  to  be  extracted  from  the 
(user-specified)  larger-scale  model. 

The  L’SER_0NLINE_TERM1NAL  interface  will  obtain  the  required 
information  via  a  standard  INPUT  and  OUTPUT  connected  "remote" 
terminal  and  the  prompting  program  will  write  that  information  to 
a  file  in  Fortran  200  NAMELIST  format  (as  used  on  the  Cyber  205, 
see  ref.  n)  for  later  transfer  and  input  to  the  appropriate  larg¬ 
er-scale  atmospheric  prediction  system. 

This  interface  will  use  free-form  query  and  response  to 
obtain  the  USER_ONLINE_TERMINAL  data  elements  presented  in  Appen¬ 
dix  D. 

3.3.2  NAMELIST_DATA.  This  interface  consists  of  the  NAMELIST 
input  to  the  NABL  Model.  (NOTE;  NAMELIST  input/output  is  a  CDC 
extension  to  ANS I -S t andard  FORTRAN  and  is  approved  for  use  at 
FNOC .  ) 

The  NAMEL I ST_DATA  interface  provides  all  of  the  unique 
details  for  a  particular  execution  of  the  NABL  model.  The  NAME- 
LIST  sequential  data  file  is  prepared  off-line.  It  will  be 
appended  to  the  batch-mode  job  which  causes  execution  of  the  NABL 
forecast  model. 

The  data  elements  contained  in  the  NAMELI  ST_DATA  file  are 
presented  in  Appendix  D. 

3.3.3  INITIAL_COND_DATA.  This  interface  provides  incoming  ini¬ 
tial  conditions  and  hourly  forecast  tendency  data  from  the  larg¬ 
er-scale  forecast  model. 
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The  IN ITIAL_COND_DATA  interface  provides  access  to  NOGAPS  or 
NORAPS  initial  condition  data  and  forecast  data  for  all  of  the 
points  specified  by  the  user  through  the  USER_ONLINE_TERMINAL 
interface  (paragraph  3.3.1)  and  passed  to  the  larger-scale  fore¬ 
cast  system  via  the  DATAPO I NT_DETA I L S  interface  (paragraph 

3.3.4) .  The  initial  condition  data  resides  on  a  formatted  se¬ 
quential  file  written  by  the  larger-scale  model  during  its  fore¬ 
cast  cycle. 

The  data  elements  contained  in  the  INITIAL_COND_DATA  file 
are  presented  in  Appendix  D. 

3.3.4  DATAPOINT_DETAILS/DATAPOINT_NAMELIST.  This  closely  relat¬ 
ed  pair  of  interfaces  provide  and  maintain  the  data  point  details 
(input  by  the  NABLS  user  through  the  USER_ONLINE_TERMINAL  inter¬ 
face  (paragraph  3.3.1))  which  must  be  read  by  the  larger-scale 
forecast  model  (NOGAPS  or  NORAPS)  in  order  to  prepare  the  INI- 
T I AL_COND_DATA  (paragraph  3.3.3)  input  to  the  NABL  model. 

The  DATAPO I NT_NAMEL 1ST  file  is  created  on  FED  and  "moved"  to 
SAM.  It  should  have  a  unique  permanent  file  name  which  identifies 
it  with  a  particular  large-scale  model  domain.  It  is  a  formatted 
sequential  file  in  FORTRAN  200  NAMELIST  format  (see  ref.  n). 

The  DATAPOINT_DETAI LS  file  will  reside  on  the  front-end  comput¬ 
er  (FED)  as  a  permanent  sequential  unformatted  file.  This  file 
serves  as  the  "master"  which  is  read  and  updated  (rewritten  and 
CATALOGed)  whenever  a  change  to  a  DATAPOINT_NAMELIST  file  is  re- 
qui red. 

The  master  DATAPO I NT_DETAILS  files  on  the  front-end  and  the 
DATAPOINT_NAMELIST  input  file  to  NOGAPS  on  SAM  will  contain  only 
the  data  elements  presented  in  Appendix  D.  For  complicated 
reasons  related  to  the  developmental  history  of  NORAPS  and  NABL, 
the  DATAPOINT_NAMELI ST  files  associated  with  NORAPS  must  have 
appended  at  their  beginning  a  set  of  20  integer,  logical,  or 
character  variables  which  control  global  (programmatic)  aspects 
of  NORAPS  model  execution.  These  variables,  which  are  unrelated 
to  NABL  and  can  vary  from  domain  to  domain,  are  not  listed  in 
the  Appendix  but  are  accounted  for  in  program  NRPPNT  (paragraph 

3.2.5)  which  now  generates  the  NAMELIST  file. 

3.3.5  REF_EFFECTS_DATA .  This  interface  provides  a  character 
(text)  display  of  refractive  effects  for  each  point  for  each  of 
the  user-specified  NABL  forecast  times  (elements  IPT  and  ITAU  in 
the  NAMELIST_DATA  file). 

The  REF_EFFECTS_DATA  interface  provides  sequential  formatted 
file  output  suitable  for  direct  input  to  any  application  program 
which  may  have  use  for  such  information  (e.g.  IREPS).  Or,  the 
file  may  be  directly  printed  or  used  as  message  text. 

This  file  is  only  generated  when  NAMEL I ST_DATA  logical 
variable  FNREFC  is  .TRUE..  It  provides  data  at  each  point  for 
each  tau  for  all  40  of  the  NABL  forecast  levels. 

On  the  following  page  is  an  example  of  an  output  file  for 
one  tau  and  point  location. 
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NABL  12  HR  FORECAST  VT  91082112Z  AT  LOCATION  32. 4N  1 2 1 . 8W 


INITIALIZED  BY 

FNOC  NOGAPS  FIELDS 

1 NCLUDES 

REFRACTION 

40  LEVELS 

SFC  M-VAL= 

362.9 

TOPO  = 

32.6 

WATER  SFC 

PRES 

TEMP 

DEWPT 

M 

GRAD  M 

HT 

REFRACTION 

(C) 

DPRSN 

1000  FT 

(FT) 

1016.8 

17.  1 

3.  1 

344.2 

49.0 

EVDT=  39FT 

1016.2 

17.1 

3.3 

343.7 

-28.905 

65.4 

TRAPP  1 NG 

1015.0 

16.9 

3.5 

343.9 

8.128 

98.2 

SUPER 

1012.6 

16.7 

3.5 

345.6 

24.924 

163.9 

NORMAL 

1009.7 

16.4 

3.5 

348  3 

32.607 

245.9 

NORMAL 

1005.5 

16.1 

3.3 

352.4 

35.844 

360.7 

NORMAL 

- 

etc.  on  up 

to  levels 

39  and 

40  - 

607.4 

2.9 

26.4 

859.3 

42.395 

14304 . 3 

NORMAL 

596.2 

1  .  7 

25.4 

880.4 

42 . 806 

14796.4 

NORMAL 

Note: 

"TOPO"  in  the 

third  line  is  terrain  height  (ft). 

Also 

,  "WATER  SFC"  would 

be  "LAND 

SFC"  if 

the  point 

was  considered  "dry"  (if  its  ground  wetness  value 
was  less  than  1.0). 


3.3.G  SKEWT_DATA.  This  interface  provides  a  character  (text) 
display  of  the  standard  "vertical  sounding"  values  required  to 
prepare  a  "SKEWT"  diagram  based  at  each  point  for  each  of  the 
user-specified  NABL  forecast  times  (elements  IPT  and  ITAU  in  the 
NAMELIST_DATA  file). 

The  SKEWT_DATA  interface  provides  standard  FNOC  sequential 
formatted  file  output  in  "English"  units  suitable  for  direct 
input  to  any  application  program  which  may  have  use  for  such 
information  (e.g.  IREPS).  Or,  the  file  may  be  directly  printed 
or  used  as  message  text. 

This  interface  also  provides  a  second  slightly  different 
"SKEWT"  file  in  metric  units  which  has  been  used  DIAGNOSTICALLY 
by  NOARL  for  input  to  a  "SKEWT"  plotting  (graphics)  program. 

The  first  file  is  only  generated  when  NA.MEL  I  ST_DATA  logical 
variable  FNSKEW  is  .TRUE,  and  the  second  when  NAMELIST  variable 
NBSKEW  is  .TRUE.  Both  provide  data  at  each  point  for  each  tau 
for  all  40  of  the  NABL  forecast  levels. 

On  the  following  page  are  examples  of  the  two  types  of 
output  -  the  same  one  tau  and  point  location  for  each  type  from 
the  same  forecast  event  -  FNOC  format  first. 
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NABL  12  HR  FORECAST  VT  91082112Z  AT  LOCATION  32. 4N  121. 8W 
INITIALIZED  BY  FNOC  NOGAPS  FIELDS  INCLUDES  WINDS 

40  LEVELS  WATER  SFC 


PRES 

TEMP 

DEWPT 

DD 

FF 

HT 

(C) 

DPRSN 

(FT) 

1016.8 

17.1 

3. 1 

325 

16 

49.0 

1016.2 

17.  1 

3.3 

325 

16 

65.4 

1015.0 

16.9 

3.5 

325 

17 

98.2 

1012.6 

16.7 

3.5 

325 

17 

163.9 

1009.7 

16.4 

3.5 

325 

18 

245.9 

1005.5 

16. 1 

3.3 

326 

18 

360.7 

etc 

.  on  up 

to  levels 

39  and 

40  - 

607.4 

2.9 

26.4 

180 

5 

14304.3 

596.2 

1  .7 

25.4 

180 

5 

14796 . 4 

Note:  Column  one  (pressure)  is  in  millibars,  wind  direction 
(DD)  is  in  degrees  true  and  wind  speed  (FF)  is  in  knots.  In  the 
third  line  "WATER  SFC"  would  read  "LAND  SFC"  if  the  point  was 
considered  "dry"  (if  its  ground  wetness  value  was  less  than  1.0). 

Also  note:  It  is  recommended  that  "TERRAIN  HEIGHT  xxx  FT"  be 
added  to  the  third  line  as  potentially  useful  information. 


2 


12789  91082 

100  12 

32.4  -121 . 

8  110 

0  1  2  119.0  82.0  1.0000 

♦  line  two  above 

cont  ’  d 

♦  9.4  .0 

17.7 

1016.78 

17.132 

3.144 

325.0 

8.53 

14 

1016.18 

17.052 

3.335 

325.0 

8.73 

19 

1014.99 

16.926 

3.472 

325.1 

8.93 

29 

1012.62 

16.702 

3.522 

325. 1 

9.16 

49 

1009.66 

16.442 

3.454 

325.3 

9.38 

74 

1005 . 52 

16.089 

3.282 

325.6 

9 . 64 

109 

etc.  on  up  to  levels  39 

and  40  - 

607.61 

2.867 

26.449 

179.7 

2.72 

4359 

596.23 

1.676 

25.448 

179.7 

2.72 

4509 

Note: 

Co  1 umns 

above  are 

pressure  (mb . ) 

,  temperature 

(C)  , 

dewpoint  depression 

(C) ,  wind 

direction  (deg. 

true) ,  wind 

speed 

(mps),  and  NABL  grid  point  height  (m)  .  Line  one  is  the  number 
of  forecast  times  requested.  Line  two  is  block-station  identifi¬ 
er,  forecast  base-DTG,  forecast  tau,  latitude,  longitude,  five 
single-integer  plot  constants,  model  source  indicator  ("2"  means 
NOGAPS),  "I"  coordinate  of  the  data  point,  "J"  coordinate, 
ground-wetness  value  (1.0000),  terrain  height  (m),  a  land-or-sea 
flag  (0  indicates  a  water  point  -  see  I N 1 T 1  AL_COND_DATA  data 
element  "PL_LRS"  in  Appendix  D  and  read  Appendix  E  for  a  complete 
description  of  this  flag),  and,  last,  surface  temperature  (C). 
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Also  Note:  A  file  nearly  identical  to  the  second  (metric) 
file  on  the  previous  page  is  produced  when  NAMEL I ST_DATA  variable 
NRSKEW  is  ".TRUE.''.  This  third  file  contains  SKEWT  data  based 
strictly  on  the  IN  I TI AL_COND_DATA  from  the  large-scale  model.  It 
permits  comparison  with  the  NABL  model  output  and  is  strictly 
DIAGNOSTIC. 

3.3.7  COMPLETE_FORECAST_DATA.  The  COMPLETE_FORECAST_DATA  inter¬ 
face  will  provide  any  (future)  interfacing  program  a  source  file 
for  any  desired  NABL  forecast  variable  -  assuming  the  resultant 
file  is  cataloged  and  made  available  for  the  appropriate  applica- 
t i ons . 

The  data  elements  to  be  contained  in  the  COMPLETE_  FORECAST_ 
DATA  file  are  presented  in  Appendix  D. 

Note:  this  interface  is  not  provided  for  in  the  existing 
NABL  code.  It  should  consist  of  a  compact  sequential  formatted 
file  containing  all  of  the  forecast  variables  for  each  forecast 
point  and  tau.  (There  is  provision  within  the  existing  code  for 
extensive,  hut  not  compact,  writes  of  many  of  the  variables  to 
standard  OUTPUT.  The  existing  OUTPUT  is  considered  strictly 
DIAGNOSTIC. ) 
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4.  DETAILED  DESIGN.  This  section  describes  the  detailed  design 
of  each  CSC.  This  is  done  by  presenting  and  discussing  each  CSU 
(program  module)  associated  with  a  particular  CSC. 

All  procedure  modules  are  written  in  CCL.  A  complete  list¬ 
ing  of  each  procedure  is  provided  in  Appendix  B. 

All  program  modules  are  written  in  AN'S  I -S  tandard  FORTRAN 
using  only  FNOC-al lowed  non-ANSI  features  (the  few  exceptions  to 
this  rule  are  noted).  A  listing  of  the  prologue  for  each  program 
and  subroutine  is  provided  in  Appendix  C.  Subroutine  modules  are 
presented  alphabetically  following  their  respective  program 
modules  -  the  same  order  as  is  used  to  present  their  prologues  in 
Appendix  C. 

Hierarchy  diagrams  are  provided  for  all  programs  and  the 
associated  calls  are  mentioned  in  the  discussion.  Calls  to 
standard  FNOC  library  subprograms,  CDC  system  calls  and  uses  of 
standard  intrinsic  functions  are  not  documented.  For  such  de¬ 
tails  the  reader  would  be  better  served  by  the  standard  CDC 
Compiler  Listings  and  Reference  Maps. 

4.1  Program  NABL  (NABL).  Figure  4-1  shows  the  hierarchy  of  the 
twenty  Program  NABL  modules. 

4.1.1  Module  NABL  (Mod_NABL).  This  is  the  main  driver  (PROGRAM) 
"module"  of  the  NABL  "component"  of  the  NABL  "system". 

Mod_NABL  is  larger  and  directly  serves  several  more  func¬ 
tions  than  would  the  normal  "driver"  in  a  more  modularized  pro¬ 
gram.  In  addition  to  initializing  control  variables  and  data 
arrays  (based  largely  on  PARAMETER  and  DATA  statements  within 
this  main  program  and  on  information  read  from  file 
NAMELIST_DATA) ,  this  module  reads  in  all  of  the  I N I T I AL_COND_ 
DATA.  The  primary  prediction  and  output  cycles  (taus  within 
geographic  points)  are  controlled  from  within  this  program  based 
on  the  NAMELIST_DATA. 

Program  NABL  calls  13  level-1  subroutines  as  shown  in  figure 

4-1  . 

Note:  It  is  suggested  that  the  reading  of  IN ITI AL_COND_DATA 
(TAPE77)  be  placed  into  a  new,  separate  subroutine.  This  could 
he  done  at  the  same  time  as  provisions  are  made  for  the  writing 
of  a  COMPLETE_FORECAST_DATA  file  as  suggested  in  paragraph  3.3.7. 
Also  at  that  time,  the  main  po i n t - by - po i n t  loop  and  the 
time-step-by-time-step  computational  loop  within  that  first  loop 
could  both  be  moved  out  into  a  separate  subr out i ne ( s ) .  These 
changes  would  add  greatly  to  the  modularity  and  therefore,  pre¬ 
sumably,  to  the  maintainability  of  this  complex  code. 

4.1.2  Subroutine  BOUND  (BOUND).  This  subroutine  is  called  by 
program  NABL.  It  computes  the  stability-dependent  surface  fluxes 
of  heat,  moisture  and  momentum;  computes  the  1  owe r -boundary 
turbulence  variables;  and  calculates  the  evaporative  duct  height. 

4.1.3  Subroutine  CORCR  (CORCR).  This  level-2  subroutine  is 
called  by  subroutine  CORRD.  It  calculates  longwave  upward  and 
downward  fluxes  and  heating  rates  and  comparable  values  for 
shortwave  (solar)  radiation. 

Subroutine  CORCR  calls  subroutines  CORWL  and  CORWS. 
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Figure  4-1.  Progran  NftBL  hierarchy. 
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4.1.4  Subroutine  CORRD  (CORRD).  This  subroutine  is  called  by 
program  NADL.  It  is  the  main  radiation  subprogram.  Integrated 
column  abundances  of  water  vapor  and  carbon  dioxide  are  computed, 
the  black  body  functions  at  each  forecast  level  are  calculated 
and  solar  flux  at  the  top  of  the  atmosphere  is  defined.  Then 
fluxes  are  calculated  by  calling  subroutine  CORCR. 

4.1.5  Subroutine  CORWL  (CORWL).  This  level-3  subroutine  is 
called  by  subroutine  CORCR.  It  computes  the  transmission  func¬ 
tions  which  are  needed  to  calculate  1 ong-wave  radiative  fluxes. 

4.1.6  Subroutine  CORWS  (CORWS).  This  level-3  subroutine  is 
called  by  subroutine  CORCR.  It  computes  the  transmission  func- 
lions  which  are  needed  to  calculate  short-wave  radiative  fluxes. 

4.1.7  Subroutine  ESAT  (ESAT).  This  subroutine  is  called  by  five 
different  level  1  subroutines  (see  figure  4-1).  It  computes  the 
saturation  vapor  pressure  (mb.)  for  a  given  temperature  (K) . 

4.1.8  Subroutine  HOSKEEP  (HOSKEEP).  This  subroutine  is  called 
by  program  NABL.  It  calculates  the  grid  volume  liquid  water 
content,  the  fractional  cloudiness,  and  the  eddy  transfer  coeffi¬ 
cients  . 

Subroutine  HOSKEEP  calls  subroutine  ESAT. 

4.1.9  Subroutine  INITIAL  (INITIAL).  This  subroutine  is  called 
by  program  NABL.  It  interpolates  the  initial  condition  data  from 
the  NOGAPS  or  NORAPS  sigma-levels  to  the  f i xed-e 1 eva t i on  points 
on  the  NABL  grid.  Also,  several  arrays  used  to  store  turbulence 
variables  are  initialized  to  zero  and  initial  values  for 
radiative  flux  transfers  above  the  computational  grid  are 
calculated. 

Subroutine  INITIAL  calls  subroutines  ESAT  and  VINTRP. 

4.1.10  Subroutine  NABL42  (NABL42).  This  subroutine  is  called  by 
program  NABL  whenever  NAMEL I ST_DATA  logical  variable  FNSKEW 
is  .TRUE..  It  writes  a  SKEWT^DATA  point-tau  logical  record  to 
TAPE42  in  FNOC  standard  (English  units)  format. 

Subroutine  NABL42  calls  subroutine  UVSD. 

4.1.11  Subroutine  NABL46  (NABL46).  This  subroutine  is  called  by 
program  NABL  whenever  NAMEL I ST_DATA  logical  variable  FNREFC 
is  .TRUE..  It  writes  a  REF_EFFECTS_DATA  point-tau  logical  record 
to  TAPE46  in  FNOC  standard  (English  units)  format. 

4.1.12  Subroutine  NABL75  (NABL75).  This  subroutine  is  called  by 
program  NABL  whenever  NAMEL I ST_DATA  logical  variable  NBSKEW 
is  .TRUE..  It  writes  a  SKEWT_DATA  point-tau  logical  record  to 
TAPE75  in  NOARL  (metric  units)  format. 

Subroutine  NABL 7 5  calls  subroutine  UVSD. 

4.1.13  Subroutine  NABL80  (NABL80).  This  subroutine  is  called  by- 
program  NABL  whenever  NAMEL I ST_DATA  logical  variable  NRSKEW 
is  .TRUE..  It  writes  a  SKEWT-type  point-tau  logical  record  to 
TAPE80  in  NOARL  (metric  units)  format.  [Note:  these  records  re- 
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fleet  the  larger-scale  model’s  output  (as  input  to  N’ABL)  and  not 
the  results  of  the  NABL  model's  forecast.] 

Subroutine  N’ABLSO  calls  subroutines  ESAT  and  UVSD. 

4.1.14  Subroutine  OUT  (OUT).  This  subroutine  is  called  by 
program  NABL.  It  writes  the  contents  of  nearly  all  forecast 
variables  to  standard  OUTPUT  (unless  NAMELI ST_DATA  logical  varia¬ 
ble  DSKIP  is  .TRUE.,  in  which  case  this  output  is  substantially 
abbreviated).  Note:  this  subroutine  could  be  replaced  or  largely 
rewritten  to  provide  for  the  COMPLETE_FORECAST_DATA  file  output 
discussed  in  paragraph  3.3.7. 

Subroutine  OUT  calls  subroutine  ESAT. 

4.1.15  Subroutine  TIIDAG  (TMDAG).  This  subroutine  is  called  by 
program  NABL  only  if  logical  NAMEL I ST_DATA  variable  TSER 
is  .TRUE..  This  routine  loads  selected  NABL  forecast  variables 
into  the  large  real  array  PDATA.  Program  NABL  writes  the  PDATA 
array  to  TAPE70  for  DIAGNOSTIC  off-line  "time  series  plotting"  at 
the  very  end  of  program  execution  (but  only  if  TSER  is  .TRUE.). 

Subroutine  TMDAG  calls  subroutine  UVSD. 

4.1.16  Subroutine  TQTEND  (TQTEND).  This  subroutine  is  called  by 
program  NABL.  It  computes  the  hourly  average  adiabatic  tenden¬ 
cies  for  temperature  and  moisture  based  on  the  1 N I T 1 AL_COND_ 
DATA. 

Subroutine  TQTEND  calls  subroutines  ESAT  and  VINTRP. 

4.1.17  Subroutine  UTENDS  (UTENDS).  This  subroutine  is  called  by 
program  NABL.  It  computes  the  hourly  average  wind  component 
tendencies  based  on  the  I N I T I AL_COND_DATA .  The  larger-scale 
model  wind  data  are  interpolated  to  the  NABL  grid  points  in  this 
subprogram. 

Subroutine  UTENDS  calls  subroutine  VINTRP. 

4.1.18  Subroutine  UVCOMP  (UVCCTilP).  This  subroutine  is  called  by 
program  NABL.  It  computes  the  momentum  flux  divergence  for  a 
given  wind  component. 

4.1.19  Subroutine  UVSD  (UVSD).  This  subroutine  is  called  by 
subroutines  NABL42,  NABL75,  NABL80  and  TMDAG.  It  converts  U  and 
V  wind  components  to  speed  (m/s)  and  direction  (deg.T). 

4.1.20  Subroutine  VINTRP  (VINTRP).  This  level-2  subroutine  is 
call  by  subroutines  INITIAL,  TQTEND  and  UTENDS.  It  interpolates 
values  from  the  larger-scale  NOGAPS  or  NORAPS  sigma-levels  to 
NABL  ver t i cal -gr id  coordinates. 

4.2  Procedure  NGPOINT  (NGPOINT).  This  procedure  executes  inter¬ 
active  program  NGNPNT  and,  if  upon  completion  of  that  program 
sense-switch  one  is  set  ON  (indicating  that  a  new 
DATAPOINTS_NAMELI ST  file  for  NOGAPS  use  has  been  created),  this 
procedure  then  executes  (BEGINs)  companion  procedure  NGNMOD. 

This  procedure  is  listed  in  Appendix  B. 

Note;  It  is  suggested  that  this  procedure  be  combined  with 
procedure  NRPNABL  below. 
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4.3  Procedure  NRPNABL  (NRPNABL).  This  procedure  executes  inter¬ 
active  program  NRPPNT  and,  if  upon  completion  of  that  program 
sense-switch  one  is  set  ON  (indicating  that  a  new 
DATAPOINTS_NAMELIST  file  for  NORAPS  use  has  been  created),  this 
procedure  then  executes  (BEGINs)  intermediary  procedure  PRO. 

This  procedure  and  procedure  PRO  are  listed  in  Appendix  B. 

Note:  It  is  suggested  that  this  procedure  be  combined  with 
procedure  NGPOINT  above,  using,  for  example,  sense-switch  two  to 
indicate  which  follow-on  transfer  procedure  (NGNMOD  or  PRO) 
should  be  executed. 

4.4  Program  NGNPNT  (NGNPNT).  Figure  4-2  shows  the  hierarchy  of 
the  nine  Program  NGNPNT  modules. 

Note:  It  is  suggested  that  this  program  be  combined  with 
program  NRPPNT  below,  this  would  eliminate  subroutine  redundancy 
in  the  FNOC  source  library,  simplify  program  maintenance,  and 
ensure  a  common  user  interface  (which  should  include  prompting 
for  land-sea  information  for  NORAPS).  (See  related  note  in  para¬ 
graph  3.2.5  and  Appendix  E.) 

4.4.1  Module  NGNPNT  (Mod_NGNPNT )  .  This  is  the  main  driver 
(PROGRAM)  "module"  of  the  NGNPNT  "component"  of  the  NABL 
"system".  In  good  modular  fashion  it  does  no  more  than  connect 
INPUT  and  OUTPUT  to  the  user’s  on-line  terminal  and  then  call 
level-1  subroutines  to  do  the  major  work.  Before  exiting,  this 
program  sets  sense  switch  one  ON  if  a  new  DATAPOINT_NAMELIST  file 
has  been  created. 

Program  NGNPNT  calls  five  level-1  subroutines  as  shown  in 
f igur e  4-2 . 

4.4.2  Subroutine  NGEXTR  (NGEXTR).  This  level-3  subroutine  is 
called  by  subroutine  NGLL.  It  extracts  characters,  integers  or 
real  numbers  from  character  strings. 

4.4.3  Subroutine  NGFORM  (NGFORM).  This  subroutine  is  called  by 
program  NGNPNT.  It  writes  the  user's  point  information  in  CDC 
FORTRAN  200  NAMELIST  format  to  TAPES  ( DATAPOI NT_NAMEL I  ST  file). 

4.4.4  Subroutine  NGLL  (NGLL).  This  level-2  subroutine  is  called 
by  subroutine  NGRPT.  It  provides  an  interface  to  the  on-line 
user  for  the  purpose  of  obtaining  valid  point  latitudes  and 
longitudes  and  prompts  for  the  user's  land  or  ocean  assessment  of 
the  point . 

Subroutine  NGLL  calls  subroutines  NGEXTR,  NGSTOP  and  NGYN. 

4.4.5  Subroutine  NGOLD  (NGOLD).  This  subroutine  is  called  by 
program  NGNPNT.  At  the  top  this  subroutine  reads  existing  (old) 
point  information  from  TAPE6  (file  NGNREC)  into  internal  arrays. 
At  alternate  entry  point  NGWRT  it  writes  the,  presumably  new  or 
revised,  point  information  back  out  to  a  new  TAPE6  (NGNREC  file). 

Subroutine  NGOLD  calls  subroutine  NGSTOP. 

Note:  Entry  point  NGWRT  is  depicted  on  figure  4-2  as  a 
separate  subroutine.  It  is  recommended  that  it  actually  be  made 
so . 
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Figure  1-2.  Progran  NGNPNT  hierarchy. 
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4.4.6  Subroutine  NGRPT  (NGRPT).  This  subroutine  is  called  by 
program  NGNPXT.  It  conducts  the  query-response  dialog  with  the 
user  -  prompting  that  person  to  enter  or  reenter  point  informa¬ 
tion  details. 

Subroutine  NGRPT  calls  subroutines  NGLL,  NGSTOP  and  NGYN. 

4.4.7  Subroutine  NGSTOP  (NGSTOP).  This  level-1,2,3  subroutine 
provides  a  common  program  exit  for  both  normal  and  abnormal 
terminal i ons . 

4.4.8  Subroutine  NGWRT  (NGWRT).  This  is  currently  an  entry 
point  in  subroutine  NGOLD.  See  paragraph  4.4.5  above. 

This  entry  calls  subroutine  NGSTOP. 

4.4.9  Subroutine  NGYN  (NGYN).  This  level-2,3  subroutine  prompts 
the  user  for  a  yes  or  no  response. 

Note:  This  subroutine  uses  "alternate  returns"  for  yes  or 

no  responses.  This  technique  is  not  permitted  by  FNOC  program¬ 
ming  standards  (ref.  k).  It  is  recommended  a  logical  variable 
return  or  something  similar  be  employed  when/if  program  NGNPNT  is 
revised.  Such  a  change  would  require  minor  changes  to  this  sub¬ 
routine  and  the  two  subroutines  which  call  it  (NGLL  and  NGRPT). 

4.5  Program  NRPPNT  (NRPPNT).  Figure  4-3  shows  the  hierarchy  of 
the  ten  Program  NRPPNT  modules. 

Note:  It  is  suggested  that  this  program  be  combined  with 

program  NGNPNT  above,  this  would  eliminate  subroutine  redundancy 
in  the  FNOC  source  library,  simplify  program  maintenance,  and 
ensure  a  common  user  interface  (which  should  include  prompting 
for  land-sea  information  for  NORAPS).  (See  related  note  in  para¬ 
graph  3.2.5  and  Appendix  E.) 

4.5.1  Module  NRPPNT  (Mod_NRPPNT) .  This  is  the  main  driver  (PRO¬ 
GRAM)  "module"  of  the  NRPPNT  "component"  of  the  NABL  "system". 
In  good  modular  fashion  it  does  no  more  than  connect  INPUT  and 
OUTPUT  to  the  user's  on-line  terminal  and  then  call  level-1 
subroutines  to  do  the  major  work.  Before  exiting  this  program 
sets  sense  switch  one  ON  if  a  new  DATAPO I NT_N AMEL I  ST  file  has 
been  created. 

Program  NGNPNT  calls  six  level-1  subroutines  as  shown  in 
figure  4-3. 

4.5.2  Subroutine  NREXTR  (NREXTR).  This  level-3  subroutine  is 
called  by  subroutine  NRLL.  It  extracts  characters,  integers  or 
real  numbers  from  character  strings. 

4.5.3  Subroutine  NRFORM  (NRFORM).  This  subroutine  is  called  by 
program  NRPPNT.  It  writes  the  user's  point  information  in  CDC 
FORTRAN  200  NAMELIST  format  to  TAPES  (the  DATAPOl NT_NAMEL I  ST 
file). 

As  previously  discussed,  this  subroutine  appends  at  the 
beginning  of  the  NAMELIST  a  set  of  20  variables  which  are  unre¬ 
lated  to  NABL.  These  variables  are  required  to  control  NORAPS 
model  execution,  and  may  vary  from  one  forecast  area  (NORAPS 
domain)  to  another.  The  values  for  these  variables  (stored  in 
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Figure  1-3.  Prograr  NRPPNT  hierarchy. 
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local  arrays  initialized  by  DATA  statements)  must  be  obtained 
from  the  person(s)  responsible  for  NORAPS  system  maintenance. 

4.5.4  Subroutine  NRLL  (NRLL).  This  level-2  subroutine  is  called 
by  subroutine  NRRPT.  It  provides  an  interface  to  the  on-line 
user  for  the  purpose  of  obtaining  valid  point  latitudes  and 
longi tudes . 

Subroutine  NRLL  calls  subroutines  NREXTR,  NRSTOP  and  NRYN. 

4.5.5  Subroutine  NROLD  (NROLD).  This  subroutine  is  called  by 
program  NRPPNT.  At  the  top  this  subroutine  reads  existing  (old) 
point  information  from  TAPE6  (file  NRNREC)  into  internal  arrays. 
At  alternate  entry  point  NRWRT  it  writes  the,  presumably  new  or 
revised,  point  information  back  out  to  a  new  TAPES  (NRNREC  file). 

Subroutine  NROLD  calls  subroutine  NRSTOP. 

Note:  Entry  point  NRWRT  is  depicted  on  figure  4-3  as  a 

separate  subroutine.  It  is  recommended  it  actually  be  made  so. 

4.5.6  Subroutine  NRPROC  (NRPROC).  This  subroutine  is  called  by 
program  NRPPNT.  It  writes  out  the  short  intermediary  procedure 
PROC  which  is  used  to  pass  the  NORAPS  area  identifier  on  into 
procedure  NRNMOD.  Procedure  PRO  is  shown  in  Appendix  B. 

Note:  NRPROC  is  the  only  subroutine  in  program  NRPPNT  which 

does  not  have  a  nearly  exact  counterpart  in  program  NGNPNT. 

4.5.7  Subroutine  NRRPT  (NRRPT).  This  subroutine  is  called  by 
program  NRPPNT.  It  conducts  the  query-response  dialog  with  the 
user  -  prompting  that  person  to  enter  or  reenter  point  informa¬ 
tion  details. 

Subroutine  NRRPT  calls  subroutines  NRLL,  NRSTOP  and  NRYN. 

4.5.8  Subroutine  NRSTOP  (NRSTOP).  This  level-1,2,3  subroutine 
provides  a  common  program  exit  for  both  normal  and  abnormal 
terminat ions . 

4.5.9  Subroutine  NRWRT  (NRWRT).  This  is  currently  an  entry 
point  in  subroutine  NROLD.  See  paragraph  4.5.5  above. 

This  entry  calls  subroutine  NRSTOP. 

4.5.10  Subroutine  NRYN  (NRYN).  This  level-2,3  subroutine 
prompts  the  user  for  a  yes  or  no  response. 

Note:  Like  its  counterpart,  NGYN ,  this  subroutine  uses 

"alternate  returns"  for  yes  or  no  responses.  This  technique  is 
not  permitted  by  FNOC  programming  standards  (ref.  k).  It  is 
recommended  a  logical  variable  return  or  something  similar  be 
employed  when/if  program  NRPPNT  is  revised  and/or  combined  with 
program  NGNPNT.  Such  a  change  would  require  minor  changes  to 
this  subroutine  and  the  two  subroutines  which  call  it  (NRLL  and 
NRRPT )  . 

4.6  Procedure  NGNMOD  (NGNMOD).  This  procedure  transfers  a 

list  of  points  ( DATAPO I NT_NAMEL I  ST )  to  SAM  and  replaces  file 
DGFILE  on  SAM  with  the  new  file.  The  job  routed  to  SAM's  input 
queue  by  this  procedure  is  OKNMD.  That  job  carries  with  it  into 
SAM  a  copy  of  the  DATAP0INT_NAMEL1 ST  file.  The  SAM  job  runs  with 


24 


NABLS  SDD/IDD  Original 


15  Novenber  1991 


operational  priority  and  requires  less  than  1  second  CPU  time. 

This  procedure  is  listed  in  Appendix  B. 

4.7  Procedure  NRNMOD  (NRNMOD).  This  procedure  is  begun  by 
executing  intermediary  procedure  PRO.  This  procedure  transfers 
a  list  of  points  (DATAPOINT_NAMELI  ST)  to  SAM  and  replaces  file 
NRNFILn  (where  integer  n  identifies  the  NORAPS  area)  with  the  new 
file.  The  job  routed  to  SAM's  input  queue  by  this  procedure  is 
OKNRD.  That  job  carries  with  it  into  SAM  a  copy  of  the  DATA- 
POINT_NAMELI ST  file.  The  SAM  job  runs  with  operational  priority 
and  requires  less  than  1  second  CPU  time. 

This  procedure  is  listed  in  Appendix  B. 
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5,  NABLS  DATA.  This  section  is  expected  to  describe  the  global 
data  elements  within  the  NABLS.  All  such  elements  are  contained 
in  the  DATAP01NT_DETAI LS  file  and/or  the  NAMEL I ST_DATA  file. 
Both  of  these  files  are  described  in  Appendix  D. 
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6.  NABLS  Data  Files.  This  section  is  expected  to  describe  each 
of  the  data  files  of  the  \ABLS  which  are  shared  by  more  than  one 
component  or  module  of  the  system.  Those  files  would  be  the 
name  1 i s t - f ormat  NAMEL I ST_DATA  file,  the  unformatted  sequential 
DATAP0INT_DETAILS  file  and  the  formatted  IN1TIAL_C0ND_F1 LE.  All 
three  of  these  files  are  described  in  Appendix  D.  (Note:  these 
files  are  really  shared  through  the  medium  of  COMMONS  into  which 
their  contents  are  read.  See  the  NABL  prologue  in  Appendix  C.) 
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7.  REQUIREMENTS  TRACEABILITY.  The  following  table  relates 
paragraphs  in  the  NABLS  requirements  document  (ref.  a)  to  the 
corresponding/supporting  paragraphs  in  this  document. 


Requirement  Design  Specification 

(reference  a.  paragraph)  (paragraph  herein) 

3.1  NABLS  external  interface  requirements 


3.1.1  USER_0NLINE_TERM1NAL .  3.3.1 

3.1.2  NAMELIST_DATA .  3.3.2 

3.1.3  INITIAL_COND_DATA .  3.3.3 

3.1.4  DATAPOINT_DETAILS/DATAPOINT_NAMELIST .  3.3.4 

3.1.5  REF_EFFECTS_DATA .  3.3.5 

3.1.G  SKEWT_DATA .  3.3.6 

3.1.7  COMPLETE_FORECAST_DATA .  3.3.7 

3.2  NABLS  capability  requirements 

3.2.1  Interface  Larger-scale  Models 

3. 2. 1.1  Interface  Control .  4.2  &  4.3 

3. 2. 1.2  Get  User  Input .  4.4  &  4.5 

3. 2. 1.3  Transfer  Points .  4.6  &  4.7 

3.2. 1.4  Transfer  Model  Ou t pu t . . ( ex t erna 1  to  NABLS) 

3.2.2  Process  Input  Data .  4.1 

3.2.3  Compute  Forecast  Variables .  4.1 

3.2.4  Prepare  Output .  4.1 

3.3  NABLS  internal  interfaces 

3.3.1  Initiate  User  Input .  4.2  &  4.3 

3.3.2  Convey  Datapoint  Details .  4.4  &  4.5 

3.3.3  Initiate  Datapoint  Transfer .  4.2  &  4.3 

3.3.4  Convey  Output  Requirements .  4.1.1 

3.3.5  Convey  Larger-Scale  Data .  4.1.1 

3.3.6  Convey  Forecast  Data .  4.1.1  &  4.1.10-15 

3.3.7  Control  Model  Execution .  4.1.1 

3.4  NABLS  data  element  requirements  (internal). 

3.4.1  SENSE_SWITCH .  4.4.1  &  4.5.1 

3.4.2  NORAPS_PROC .  4.5.6  &  4.7 

3.4.3  DATAPOINT_NAMELIST .  4.4.3  &  4.5.3 

3.4.4  OUTPUT_CONTROL_FLAGS .  4.1.1 

3.4.5  LARGER_SCALE_DATA .  4.1.1 

3.4.6  FORECAST  VARIABLES .  4.1.1  &  4.1.10-15 
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8.  Notes. 


8.1  Acronyms,  abbreviations  and  identifiers. 


This  appendix  defines  acronyms  and  abbreviations  used  in  the 
document  text.  It  also  provides  a  long-title  and  paragraph 
reference  for  each  project  unique  identifier  (except  for  subrou¬ 
tine  modules). 


_ Short  Title 

CCL 

CDC 

COMPLETE_FORECAST_ 

DATA 

CSC 

CSCI 

CSC 

DATAP0INT_DETAILS 

D AT APO 1 NT_N  AM  E  L I S  T 

DFD 

FED 

FNOC 

IDD 

INITIAL_COND_ 

DATA 

IREPS 

IRS 

NABL 

NABLS 

NAMELIST_DATA 

NGNMOD 


Long  Title  and/or  Explanation _ 

(CDC)  Cyber  (job)  Control  Language 
Control  Data  Corporation 
Complete  Forecast  Data  External 
Interface/File  (para.  3.1.7) 

Computer  Software  Component  (subset  of  a 
CSCI  ) 

Computer  Software  Configuration  Item 
Computer  Software  Unit  (subset  of  a  CSC  ) 
Datapoint  Details  External  Interface/File 
and  Internal  Data  Set  (para.  3.3.4) 
Datapoint  Namelist  Internal  Interface/File 
( para .  3.3.4) 

Data  Flow  Diagram 

The  designated  "Front-End"  computer 
(usually  a  CDC  Cyber  860) 

Fleet  Numerical  Oceanography  Center, 
Monterey,  CA 

Interface  Design  Document 
Initial  Conditions  Data  External 
Interface/File  (para.  3.3.3) 

Integrated  Refractive  Effects  Prediction 
System 

Interface  Requirements  Specification 
Naval  Atmospheric  Boundary  Layer  (model  or 
program) 

NABL  Model  System 

Namelist  Data  External  Interface/File 
( par a .  3.1.2) 

NOGAPS  Transfer  Points  Procedure  (para. 
3.2.6  and  4.6) 
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_ Short  Title 

NGNPNT 

NGPOINT 

NOARL 

NOGAPS 

NORAPS 

NOS/BE 

NRNMOD 

NRPNABL 

NRPPNT 

REF_EFFECTS_DATA 

SAM 

SDD 

SKEWT_DATA 

SRS 

L'SER_ONLINE_ 

TERMINAL 

VSOS 


_ Long  Title  and/or  Explanation _ 

NOGAPS  Get  User  Input  Program  (para.  3.2.4 
and  4.4) 

NOGAPS  Interface  Control  Procedure  (para. 

3.2.2  and  4.2) 

Naval  Oceanographic  and  Atmospheric 
Research  Laboratory 
Navy  Operational  Global  Atmospheric 
Prediction  Model 

Navy  Operational  Regional  Atmospheric 
Prediction  Model 

(CDC’s)  Network  Operating  System/Batch 
Env i ronment 

NORAPS  Transfer  Points  Procedure  (para. 
3.2.7  and  4.7) 

NOGAPS  Interface  Control  Procedure  (para. 

3.2.3  and  4.3) 

NORAPS  Get  User  Input  Program  (para.  3.2.5 
and  4.3) 

Refractive  Effects  Data  External 
Interface/file  (para.  3.3.5) 

The  CDC  Cyber  205  "super-computer"  System 
Software  Design  Document 
SKEWT  Data  External  Interface/File 
(para.  3.3.6) 

Software  Requirements  Specification 
User  Online  Terminal  External  Interface 
(para.  3.3.1) 

(CDC's)  Virtual  Storage  Operating  System 
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Appendixes . 


Appendix  A.  Data  flow  diagram  (DFD)  symbology 

A  data  flow  diagram  (DFD)  is  a  graphical  technique  used  to  depict 
information  flow  and  the  transforms  that  are  applied  as  data  move 
from  input  to  output.  It  is  also  known  as  a  bubble  chart.  The 
following  illustrates  the  basic  form  of  a  DFD  as  utilized  in  this 
document . 


External 

Entity 


A  source  of  system  inputs,  or  sink  of  system 
outputs. 


Performs  some  transformation  of  its  input  data 
to  yield  its  output  data. 


Data  Flow  Shows  the  data  items  moving  into  or  out  of 

- >  of  processes,  data  stores,  sources  or  sinks; 

the  arrowhead  indicates  the  direction  of  data 
t  rans  f  er  . 


Control  Flow  Shows  the  exercise  of  control  of  data  flow 

———————  >  from  one  exercised  by  one  external  entity  or 

process  upon  another;  the  arrowhead  indicates 
direction  of  control  transfer. 


-  A  repository  of  data,  in  a  logical  sense;  it 

Data  Store  can  be  implemented  as  a  disk  file,  a  memory 

-  buffer,  or  a  database,  etc. 
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Appendix  B.  Procedure  File  Listings 

Procedure  Name  Page 


NGPOINT 

B-2 

NRPNABL 

B-3 

NGNMOD 

B-4 

PRO 

B-5 

NRNMOD 

B-5 

B-1 
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Listing  of  Procedure  NGPOIN’T. 


NGPOINT. 

LIBRARY(*FNWCL1B) 

APLIB( ,*KEEP) 

NDF.OFF . 

KEEP(FNWCLIB) 

ETL, 100. 

MOUNT, SNrSHARDSK . 

APLIB(MT1731 ,*NGNPNT) 

NDF.ON. 

OFFSW( 1 ) 

NGNPNT. 

RETURN . TAPE7 , NGNPNT . 

COMMENT.  TEST  FOR  AND  RUN  PRO  HERE 
I FX(SSW, 1 ,NOSAM) 

NDF ,OFF . 

REWIND, TAPES. 

COPYSBF , TAPES, OUTPUT. 

REWIND, TAPES. 

NDF , ON . 

BEGI N,NGNMOD,OPSF 1 L . 

ENDIF(NOSAM) 

DISCONT, INPUT. 

D I SCONT, OUTPUT. 

NDF, OFF. 

APL1B( ,*KEEP) 

KEEP(FNWCL IB) 

DSMOUNT , SN=SHARDSK . 

NDF, ON. 

REVERT , 

EX!T(P) 

Dl SCONT( I NPUT) 

Dl SCONT(OUTPUT) 

APL(B( ,*KEEP) 

KEEP(FNWCL IB) 

DSMOUNT , SN=SHARDSK . 

NDF , ON . 

REVERT. 

EXIT. 

REVERT(ABORT) 

EXIT. 

EXIT. 

C  INTERFACE  PROCEDURE  TO  SET  NEW  NABL  LOCATIONS  ON  SAM  FILE 
C  DGFILE  FOR  NOGAPS. 

C  LEWIT  CSC  12/90 
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Listing  of  Procedure  NRPNABL. 


NRPNABL. 

L1BRARY(*FNWCL1B) 

APLIB( ,*KEEP) 

NDF,OFF. 

KEEP(FNWCL IB) 

ETL, 100. 

MOUNT, SN=SHARDSK . 

APLIB(MT1731 .♦NRPPNT) 

NDF,ON. 

OFFSW( 1 ) 

NRPPNT. 

COMMENT.  TEST  FOR  AND  RUN  PRO  HERE 
IFX(SSW, 1 .NOSAM) 

NDF.OFF. 

REWIND, TAPES. 

COPYSBF , TAPES , OUTPUT . 

REWIND, TAPES . 

NDF,ON. 

REWIND, PRO. 

BEGIN, PRO, PRO. 

ENDIF(NOSAM) 

DISCONT, INPUT. 

D I SCONT, OUTPUT. 

NDF,OFF. 

APLI8( ,^KEEP) 

KEEP(FNWCL IB) 

DSMOUNT , SN=SHARDSK . 

NDF,ON. 

REVERT. 

EXIT(P) 

D I SCONT ( INPUT) 

D I SCONT(OUTPUT) 

APLIB( ,*KEEP) 

KEEP(PNWCL IB) 

DSMOUNT, SN=SHARDSK. 

NDF ,ON . 

REVERT . 

EXIT. 

REVERT(ABORT) 

EXIT. 

EXIT. 

C  INTERFACE  PROCEDURE  TO  SET  NEW  NABL  LOCATIONS  ON  SAM  FILE 
C  NRNFIL7  FOR  NOR  -  S. 

C  LEW  IT  CSC  12/5:';. 
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Listing  of 


Procedure  NGNMOD. 


NQNMOD. 

COPYBR , SAMJOB , SAM  I N , 1 . 

COPYBF, TAPES, SAMI N. 

RETURN, SAMJOB. 

REWIND,SAMIN. 

ROUTE, SAM  IN, DC= IN. 

REVERT. 

EXIT. 

REVERT, ABORT. 

EXIT. 

EXIT. 

PROCEDURE  TO  MAKE  NEW  VERSION  OF  NABL  INPUT  FILE  DGF I LE  ON 
SAM  FOR  NOGAPS.  INITIATED  BY  PROCEDURE  NGPOINT. 

CHANGE  RECORD 

<<CHANGE  NOTICE>>  NGNMOD*01  (  20  MARCH  1991)  --  LEW  IT 
ROUTPUT  OUTPUT  TO  BIN  46 
CREATED  BY  H.  LEWIT  CSC  12/90 
/DATA  SAMJOB 
OKNMD,STSAM.46  LEWIT 
USER(U  =  xxxxxx , AC=xxxxxxx  ) 

RESOURCE(TL=100,PR IO=12,WS=1000,LP=5, JCAT=OPS) 

ROUTPUT ( DC=PR ,T I D=C, B I N=46 , S=0 ,CO=$DGF I LE  CHG$) 

COPY, I NPUT,DGHOLD. 

COPY, DGHOLD, OUTPUT. 

PURGE, DGF I LE. 

RETURN, DGF I LE. 

SW I TCH , DGHOLD , DGF I LE . 

DEF I NE ,DGF I LE . 

PERM  IT, DGF I LE , USER  =  * , ACCESS  =  R . 

COPY, DGF I LE, OUTPUT. 

EXIT. 
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Listing  of  Procedure  PRO. 


PRO. 

BEGlN,NRNMOD,OPSFIL,n. 

REVERT. 

Note:  the  integer  "n"  in  line  two  above  identifies  a 
particular  NORAPS  area. 


Listing  of  Procedure  NRNMOD. 


NRNMOD.NAREA. 

COPYBR , SAMJOB , SAM  I N , 1 . 

COPYBF,TAPE8,SAMIN. 

RETURN, SAMJOB. 

REWIND.SAMIN, 

ROUTE , SAMI N ,DC= I N . 

REVERT. 

EXIT. 

REVERT, ABORT. 

EXIT. 

EXIT. 

PROCEDURE  TO  MAKE  NEW  VERSION  OF  NABL  INPUT  FILE 
NRNFIL_NAREA  ON 

SAM  FOR  NORAPS.  INITIATED  BY  PROCEDURE  NRPNABL . 

CHANGE  RECORD 

CREATED  BY  H.  LEWIT  CSC  12/90 
/DATA  SAMJOB 
OKNRD,STSAM.46  LEWIT 
USER(U=xxxxxx,AC=xxxxxxx) 

RESOURCE(TL=100,PR IO=12,WS=1000,LP=5, JCAT=OPS) 

COPY, INPUT, NRNF IL. 

PURGE, NRNFIL_NAREA. 

SW I TCH , NRNF I L , NRNF I L_NAREA . 

DEF I NE ,NRNF I L_NAREA. 

PERMIT,NRNFIL_NAREA,USER=*,ACCESS=R. 

COPY , NRNF I L_NAREA , OUTPUT . 

EXIT. 
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Appendix  C.  Module  Proloques 


Proeram 

Modu 1 e 

Level 

Paee 

NABL 

NABL 

0 

C-2 

BOUND 

1 

C-5 

CORCR 

2 

C-6 

CORRD 

1 

C-7 

CORWL 

3 

C-8 

CORWS 

3 

C-9 

ESAT 

2 

C-10 

HOSKEEP 

1 

C-11 

INITIAL 

1 

C-12 

NABL42 

1 

C-14 

NABL46 

1 

C-15 

NABL 7 5 

1 

C-1  6 

NABL80 

1 

C-17 

OUT 

1 

C-18 

TMDAG 

1 

C-19 

TQTEND 

1 

C-20 

UTENDS 

1 

C-21 

UVCOMP 

1 

C-22 

UVSD 

2 

C-23 

VINTRP 

2 

C-24 

NGNPNT 

NGNPNT 

0 

C-25 

NGEXTR 

3 

C-27 

NGFORM 

1 

C-29 

NGLL 

2 

C-30 

NGOLD 

1 

C-32 

NGRPT 

1 

C-33 

NGSTOP 

1,2,3 

C-34 

NGWRT 

1 

( ENTRY 

in  subroutine 

NGOLD) 

NGYN 

2,3 

C-35 

NRPPNT 

NRPPNT 

0 

C-36 

NREXTR 

3 

C-38 

NRFORM 

1 

C-40 

NRLL 

2 

C-42 

NROLD 

1 

C-43 

NRPROC 

1 

C-44 

NRRPT 

1 

C-45 

NRSTOP 

1,2,3 

C-46 

NRWRT 

1 

( ENTRY 

in  subroutine 

NROLD) 

NRYN 

2,3 

C-47 

C-1 
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PROGRAM  NABL 

. START  PROLOGUE . 

PROGRAM  NAME:  NABL 
DESCRIPTION: 

THIS  PROGRAM  INITIALIZES  A  ONE-DIhENSIONAL  BOUNDARY  LAYER 
MODEL  WITH  PROFILES  OF  WIND,  TEMPERATURE,  AND  MOISTURE 
FROM  NOQAPS  OR  NORAPS  AT  PRESELECTED  LOCATIONS.  IT  FORECASTS 
THE  EVOLUTION  OF  THESE  PROFILES  DUE  TO  PARAMETERIZED  PHYSICAL 
PROCESSES,  SURFACE  FORCING,  AND  SYNOPTIC-SCALE  ADIABATIC, 
ADVECTIVE  TENDENCIES.  HIGH-RESOLUTION  FORECASTS  OF  BOUNDARY 
LAYER  STRUCTURE  OUT  TO  24-  36  H  ARE  PRODUCED. 

ORIGINAL  PROGRAM-ER,  DATE:  STEPHEN  BURK  (NOARL)  1987 

CURRENT  PROQRATtlER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW IT  (CSC) 

LIANA  ZAMBRESKY  (FNOC) 

COMPUTER/OPERATING  SYSTEM:  CYBER  205  VSOS 

CLASSIFICATION:  UNCLASSIFIED 

INPUT  FILES: 

NAMELIST  =  INPUT  FROM  POOL  FILE  CONTAINING  DIRECTIVES  FOR 
LOCATIONS  AND  LENGTH  OF  NABL  FORECASTS  (TAPES) 

NLTAUF  -  FINAL  TAU  OF  NABL  FORECAST 

NTAUF  -  n  OF  OUTPUT  TAU(S)  FOR  NABL 

I TAU  -  TIME  (IN  HOURS)  FOR  NABL  FORECAST  OUTPUT 

NFCST  -  NUMBER  OF  LOCATIONS  FOR  NABL  FORECASTS 

IPT,JPT  -  PAIRS  OF  l,J  INPUT  GRID  POINTS 

ALAT,ALON  -  LATITUDE/LONGITUDE (0-360  W)  INPUT 

USE  1/10  DEG  ACCURACY,  AS  IN  NORAPS  INPUT 
ITYPE  -  "IJ"  (IPT,JPT)  OR  "LL"  (ALAT,ALON)  INPUT  TYPE 
NBLKST  -  5-DIQIT  BLOCKSTATION  IDENTIFIER  (OPTIONAL) 

NBSKEW  -  TRUE  IF  NABL  SKEW-T  PLOTS  DESIRED 

NRSKEW  -  TRUE  FOR  NORAPS  SKEW-T  PLOTS  (AT  12-HR  INTERVALS) 

FNSKEW  -  TRUE  FOR  FNOC  FORMAT  SKEWTS 

FNREFC  -  TRUE  FOR  FNOC  REFRACTIVITY  FILE 

NRTAU  -  TAU  AT  WHICH  NORAPS  SKEW-T  WANTED 

NSKEW  -  #  OF  NORAPS  SKEW-T  TIMES  WANTED 

DSKIP  -  TRUE  IF  LIMITED  #  OF  FIELDS  SENT  TO  OUTPUT 

TSER  -  TIME  SERIES  DIAGNOSTIC  OPTION 

TAPE77  =  INPUT  FILE  CONTAINING  INITIAL  VALUES  AND  ADIABATIC 
TEMPERATURE  AND  MOISTURE  TENDENCIES  FROM  NORAPS. 

OUTPUT  FILES: 
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C  TAPE42  =  FORECAST  DATA  FOR  SKEVfT  IN  FNOC  FORMAT 

C  TAPE46  =  FORECAST  REFRACTIVITY  VALUES 

C  TAPE75  =  NABL  FORECAST  SKEW  FILE 
C  TAPE80  =  NORAPS  FORECAST  SKEWT  FILE 
C  TAPE70  =  PDATA  -  TIME  SERIES  DIAGNOSTICS 
C 

C . MAINTENANCE  SECTION . 

C 

C  BRIEF  DESCRIPTION  OF  PROGRAM  MODULES: 

C 

C  INITIAL  -  SET  UP  INITIAL  NABL  PROFILE  FROM  NORAPS/NOQAPS  INPUT 

C  BOUND  -  SURFACE  FLUXES,  TKE  EQUATION,  EVAP  DUCT  CALCULATION 

C  HOSKEEP  -  EDDY  TRANSFER  COEFFICIENTS,  SFC  ENERGY  BUDGET 
C  -  FRACTIONAL  CLOUDINESS 

C  UVOOhf*  -  MOrtNTUM  FLUX  DIVERGENCE,  NUDGING  OF  WIND  COMPONENTS 
C  OORRD  -  MAIN  RADIATION  SUBROUTINE 

C  OORCR  -  RADIATIVE  HEATING  RATES 

C  OORWL  -  LONGWAVE  TRANSMISSION  FUNCTIONS 

C  CORWS  -  SHORTWAVE  TRANSMISSION  FUNCTIONS 

C  BLKDATA  -  INITIALIZE  CONSTANTS  AND  VARIABLES 

C  VINTRP  -  VERTICAL  INTERPOLATIONS 

C  ESAT  -  SATURATION  VAPOR  PRESSURE 

C  OUT  -  PRINT  NABL  FORECAST  OUTPUT 

C  UTENDS  -  COMPUTE  WIND  TENDENCIES  FROM  NORAPS/NOQAPS  INPUT 

C  TQTEND  -  COMPUTE  TEMPERATURE  &  MOISTURE  TENDENCES 

C  UVSD  -  CONVERT  U  &  V  W I  NO.  OC»f>ONENTS  TO  SPEED/DIRECTION 

C  NABL75  -  WRITE  NABL  FORECAST  TO  FILE  FOR  SKEW-T  PLOT 

C  NABL80  -  WRITE  NORAPS  FORECAST  TO  FILE  FOR  SKEW-T  PLOT 

C  (NOT  AVAILABLE  IN  OPERATIONAL  RUN) 

C  NABL42  -  WRITE  NABL  FORECAST  TO  FILE  IN  FNOC  FORMAT 

C  NABL46  -  WRITE  NABL  REFRACTIVITY  FORECAST  TO  OUTPUT  FILE 

C  TTDAG  -  WRITE  TO  OUTPUT  FILE  FOR  TIFE  SERIES  PLOT 
C  (FOR  DIAGNOSTIC  PURPOSES  ONLY) 

C 

C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  ABT  =  ABSOLUTE  TEMP . ( KELV I N ) 

C  ABTH  =  POTENTIAL  TEMP. 

C  DT  =  TIME  STEP (SECONDS) 

C  NN  =  TOTAL  NO.  OF  GRID  PTS. 

C  IT  =  ITERATION  NO. (COUNTER) 

C  ITOT  =  LAST  ITERATION  OF  NABL  FORECAST 

C  P(l)  =  PRESSURE(MB) 

C  QL  =  LIQUID  WATER  OONTENT(NONDI MENS lONAL.N.D. ) 

C  QW  =  TOTAL  WATER  CONTENT ( N, D. ) 

C  QWZ  =  SURFACE(SFC)  QW  VALUE 

C  SPECHUM  =  SPECIFIC  HUMIDITY 

C  TH  =  VIRTUAL  POTENTIAL  TEMP. 

C  THL  =  LIQUID  WATER  POTENTIAL  TEMP. (KELVIN) 

C  THLZ  =  SFC  VALUE  OF  THL 

C  U  =  U-COMPONENT  OF  WIND(X-DIRECT!ON) 

C  UG  =  X-OOFf>ONENT  OF  QEOSTROPHIC  WIND 

C  V  =  Y-OOMPONENT  OF  WIND 

C  VG  =  Y-OOMP.  OF  QEOSTROPHIC  WIND 

C  XMZ  =  SFC  XM  VALUE 
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C  Z(l )  =  GRID  PT.  HEIGHTS 

C  Z1  =  SFC  ROUGHNESS  HT. 

C 

C  OOrtlON  BLOCKS: 

C 

BLANK  OOWON  -  SEE  BELOW. 

ECSBLK  -  EXTENDED  CORE  MEMORY  FOR  NORAPS/NOQAPS  TENDENCIES 
BLK2  -  ARRAYS  CONTAINING  GRID  POINT  INITIAL  CONDITIONS 
AND  SPECIFICATIONS 

BLK3  -  ARRAYS  CONTAINING  HOURLY  TENDENCY  INFORMATION 
LANGUAGE:  FTN5 

RECORD  OF  CHANGES: 


END  PROLOGUE 


BLANK  OOrtlON  DESCRIPTION: 

THE  BLANK  OOTtlON  BLOCK  CONTAINS  MANY  OF  THE  MAJOR  MODEL 
VARIABLES.  THE  ARRAYS  ARE  DIMENSIONED  AOOOROING  TO  THE 
NUMBER  OF  MODEL  GRID  POINTS  IN  TTC  VERTICAL.  VARIABLE 
NAMES  ARE  USUALLY  CHOSEN  TO  BE  SUGGESTIVE  OF  THE 
MATHEMATICAL  FORM  THE  EXPRESSION  TAKES  IN  AN  EQUATION. 

FOR  EXAMPLE,  THE  SCALING  PARAFETERS  OF  MOMENTUM,  POTENTIAL 
TEMPERATURE,  AND  TOTAL  MOISTURE  ARE  WRITTEN  AS  USTAR, 

TH.STAR,  AND  QWSTAR,  RESPECTIVELY  -  SINCE  TF€SE  QUANTITIES 

ARE  WRITTEN  MATHEMATICALLY  WITH  THE  SUBSCRIPT  "STAR".  AS 
ANOTHER  EXATELE,  THE  VARIABLE  NAME  FOR  VIRTUAL  HEAT  FLUX 
IS  THVWB  BECAUSE  THIS  QUANTITY  GENERALLY  IS  WRITTEN 
MATHEMATICALLY  AS  THE  PRODUCT  OF  FLUCTUATIONS  IN  VIRTUAL 
TEMPERATURE  (THV)  AND  VERTICAL  VELOCITY  (W) ,  WITH  AN 
OVERBAR  (B)  -  HENCE,  T>C  EXPRESSION  THVWB. 

PARAMETER  (NN=40,NR=21 ,NL=1 1 ,NT=37,NP=20) 

NN  =  NUMBER  OF  NABL  VERTICAL  LEVELS 

NR  =  NUMBER  OF  NORAPS  LEVELS 

NL  =  NUMBER  OF  TENDENCY  LEVELS 

NT  =  LENGTH  OF  FORECAST  IN  HOURS  +  1  (TAU  0) 

NP  =  NUMBER  OF  OUTPUT  POINTS 
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SUBROUTINE  BOUND 


. START  PROLOGUE . 

SUBPROGRAM  NAME:  BOUND 

DESCRIPTION; 

THIS  SUBPROGRAM  COMPUTES  THE  STABILITY-DEPENDENT  SURFACE 
FLUXES  OF  HEAT,  MOISTURE,  AND  MChENTUM.  THE  ALGOR ITHiS 
DESCRIBED  BY  LOUIS,  BOUND.  LAYER  METEOR.,  17,  PP. 187-202, 

1979  ARE  USED.  ALSO,  LOWER  BOUKDARY  OOOITIONS  ON  TURBULENCE 
VARIABLES  ARE  COMPUTED  BY  ASSUMPTION  OF  A  PRODUCTION- 
DISSIPATION  EQUILIBRIUM  IN  THE  TKE  EQUATION.  EVAPORATIVE 
DUCT  HEIGHT  CALCULATIONS  ARE  ALSO  PERFORMED  HERE. 

ORIGINAL  PROGRA^t1ER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAM-ER:  ROLF  LANaAND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  2AMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  BOUND 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  OORCR(QR,FPO,FMH,TR,X1U,X2U,GAM2,EO,EH, 


1  A.PRATIO.QL, IM1,2,C1,C2,C3,WK,L2,ENU,LL, 

2  FM.FP.FN) 

C 

C . START  PROLOGUE . 

C 


C  SUBPROGRAM  NAME:  OORCR 
C 

C  DESCRIPTION: 

C 

C  THIS  SUBPROGRAM  IS  USED  IN  SEPARATE  CALLS  TO  COMPUTE, 

C  (1)  LONGWAVE  UPWARD  AND  OOWTMARO  FLUXES  AND  HEATING 
C  RATES,  AND  TVCN  (2)  THE  COMPARABLE  VALUES  FOR  SHORTWAVE 
C  (SOLAR)  RADIATION.  THE  TRANSMISSION  FUNCTION,  WHICH 
C  IS  A  FUNCTION  OF  THE  INTEGRATED  VAPOR,  ETC.,  IN  THE 
C  OOLUW,  IS  OONPUTED  BY  CALLS  TO  TVE  SUBROUTINES 

C  OORWL  AND  CORNS.  NOTE  THAT  CORWL  AND  CORNS  ARE 

C  EQUIVALENT  WITH  "TR"  THROUGH  THE  ARGUMENT  LIST.  AT  CLOUD 
C  TOP,  A  SPECIAL  ONE-SIDED  DIFFERENCE  IS  TAKEN  TO  AVOID 
C  ARTIFICIALLY  STRONG  COOLING  RATES  AT  THE  FIRST  GRID 
C  POINT  IN  THE  CLEAR  AIR  ABOVE  CLOUD  TOP. 

C 

C  ORIGINAL  PROGRAI^R,  DATE:  STEPHEN  BURK  (NOARL),  1987 

C 

C  CURRENT  PROGRA^^ER:  ROLF  LANGLAND  (NOARL) 

C  HOWARD  LEW IT  (CSC) 

C  LIANA  2AMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE): 

C 

C  CALL  OORCR(QR,FPO,FrH,TR,X1U,X2U,GAM2,EO,EH,A,PRATIO,QL, 

C  IM1 ,Z,C1,C2,C3,WK,L2,ENU,LL,FM,FP,FN) 

C 

C  INPUT  PARAMETERS: 

C 

C  FPO  =  THE  UPWARD  RADIATIVE  FLUX  FROM  THE  SURFACE 

C  FMH  =  THE  DOWNWARD  FLUX  AT  THE  "TOP"  OF  THE  ATMOSPHERE 

C  TR  =  THE  TRANSMISSION  FUNCTIOJ  COMPUTED  IN  CORWL/CORWS 
C  X1U  =  THE  "UPPER  ATMOSPHERE"  WATER  VAPOR  PATH  INTEGRAL 

C  X2U  =  THE  ATMOSPHERIC  DENSITY  INTEGRAL 

C  GAM2  =  THE  SURFACE  ALBEDO 

C  EO  =  THE  UPWARD  EMITTED  FLl«  FROM  THE  SURFACE 

C  EH  =  THE  DOWNWARD  EMITTED  FLUX  FROM  THE  TOP  OF  ATMOS. 

C  A  =  AN  ABSORPTION  COEFFICIENT  FOR  CLOUD  DROPLETS 

C  PRATIO  =  THE  EXNER  FUNCTION 

C  QL  =  THE  LIQUID  WATER  CONTENT 

C  IM1  =  NN  -  1 

C  Z  =  ARRAY  OF  GRID  POINT  HEIGHTS 

C  C1,C2,C3  =  ARRAYS  FOR  OOM»UTINQ  VERTICAL  DERIVATIVES 
C  WK  =  GENERAL  PURPOSE  WORK  ARRAY 

C  L2  =  INDEX  DETERMINES  WHETHER  WK( I ,3)  OR  WK( I ,4)  USED 

C  ENU  =  UPPER  ATMOSPHERE  BLACK  BODY  FUNCTION  VALUES 

C  LL  =  NUMBER  OF  "UPPER  ATMOSPHERE"  POINTS  —SEE  S/R  INITIAL 
C 
C 
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OUTPUT  PARAMETERS: 

OR  =  THE  COMPUTED  RADIATIVE  HEATING/  COOLING  RATES 
FM  =  DOWMARD  RADIATIVE  FLUX 

FP  =  UPWARD  RADIATIVE  FLUX 

FN  =  NET  FLUX  AT  EACH  GRID  LEVEL  (=FP( I )-FM( I ) ) 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 


SUBROUTINE  CORRD 
C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  OORRD 
C 

C  DESCRIPTION: 

C 

C  THIS  IS  THE  MAIN  RADIATION  SUBPROGRAM.  FIRST  INTEGRATED 
C  COLUMN  ABUNDANCES  OF  WATER  VAPOR  CARBON  DIOXIDE  ARE 

C  COMPUTED,  ALONG  WITH  THE  BLACK  BODY  FUNCTION  AT  EACH  LEVEL. 
C  THE  SOLAR  FLUX  AT  THE  TOP  OF  THE  ATMOSPHERE  IS  DEFINED 
C  BASED  ON  THE  LOCAL  HOUR  ANGLE  AND  SOLAR  DECLINATION.  THEN 

C  CALLS  TO  S/R  CORCR  ARE  MADE:  FIRST  TO  GET  THE  LONGWAVE 

C  FLUX  AND  HEATING  RATES,  THEN  TO  GET  THE  SHORTWAVE  VALUES. 

C  MUCH  OF  THIS  CODING  IS  BASED  ON  A.R.A.P.  REPORT  NO.  289, 

C  1976. 

C 

C  ORIGINAL  PROGRAMMER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

C 

C  CURRENT  PROGRAMMER:  ROLF  LANGLAND  (NOARL) 

C  HOWARD  LEW  IT  (CSC) 

C  LIANA  ZAMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  CORRD 
C 

C . MAINTENANCE  SECTION . 

C 

C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  OORWL(XL1 ,XL2,XL3,AL,WK,N) 
. START  PROLOGUE... 


SUBPROGRAM  NAME :  OORWL 
DESCRIPTION: 

THIS  SUBPROGRAM  COMPUTES  THE  TRANSMISSION 
FUNCTIONS  FOR  CALCULATING  RADIATIVE  FLUXES. 

OORWL  IS  FOR  LONGWAVE  TRANSMISSION  FUNCTIONS,  WHILE 
OORWS  IS  FOR  THE  SHORTWAVE  TRANSMISSION  FUNCTIONS. 

ORIGINAL  PROGRATTER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRA^^ER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  ZANBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  CORWL(XL1 ,XL2 ,XL3, AL ,WK ,N) 

INPUT  PARAMETERS: 

XL1  =  DENSITY  PATH  INTEGRAL  OF  WATER  VAPOR  FOR  TRANSMISSION  FCTN. 
XL2  =  DENSITY  INTEGRAL  OF  002  (LONGWAVE)  OR  AEROSOL  (SHORTWAVE) 
XL2  =  LIQUID  WATER  DENSITY  INTEGRAL 
AL  =  CLOUD  WATER  ABSORPTION  COEFFICIENT 
N  =  NUMBER  OF  GRID  POINTS 

OUTPUT  PARAMETERS: 

WK  =  LONGWAVE  TRANSMISSION  FUNCTION 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  OORWS  (XL1 , XL2 ,XL3, AL ,WK ,N) 
. START  PROLOGUE - 


SUBPROGRAM  NAME:  CORWS 
DESCRIPTION: 

THIS  SUBPROGRAM  COMPUTES  THE  TRANSMISSION 
FUNCTIONS  FOR  CALCULATING  RADIATIVE  FLUXES. 

CORWL  IS  FOR  LONGWAVE  TRANSMISSION  FUNCTIONS,  WHILE 
OORWS  IS  FOR  THE  SHORTWAVE  TRANSMISSION  FUNCTIONS. 

ORIGINAL  PROGRATTER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAMMER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW IT  (CSC) 

LIANA  2AMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  CORWS(XL1 ,XL2,XL3,AL,WK,N) 

INPUT  PARAMETERS: 

XL1  =  DENSITY  PATH  INTEGRAL  OF  WATER  VAPOR  FOR  TRANSMISSION  FCTN. 
XL2  =  DENSITY  INTEGRAL  OF  002  (LONGWAVE)  OR  AEROSOL  (SHORTWAVE) 
XL2  =  LIQUID  WATER  DENSITY  INTEGRAL 
AL  =  CLOUD  WATER  ABSORPTION  COEFFICIENT 
N  =  NLMBER  OF  GRID  POINTS 

OUTPUT  PARAMETERS: 

WK  =  SHORTWAVE  TRANSMISSION  FUNCTION 

. . MAI  NTENANCE  SECT  ION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  ESAT(ES1 ,T1 ,N) 

. START  PROLOGUE . 

SUBPROGRAM  NAME :  ESAT 

DESCRIPTION: 

THIS  SUBPROGRAM  OOMPUTES  THE  SATURATION  VAPOR  PRESSURE  (MB) 
FOR  A  GIVEN  TENPERATURE  (KELVIN)  BASED  ON  LOWE  (JAM,  16, 
100-103,  1977).  FOR  TEWERATURE  LESS  THAN  223K,  OONSTANTS 
ARE  USED  BASED  ON  NEPRF  REPORT  BY  LOWE. 

ORIGINAL  PROQRAltCR,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRATPER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW IT  (CSC) 

LIANA  ZAMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  ESAT(ES1,T1) 

INPUT  PARAMETERS: 

T1  =  TEMPERATURE  (KELVIN) 

N  =  NUMBER  OF  GRID  POINTS 

OUTPUT  PARAMETERS: 

ESI  =  SATURATION  VAPOR  PRESSURE  (MB) 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  SECTION 
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SUBROUTINE  HOSKEEP 


C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  HOSKEEP 
OC 

C  DESCRIPTION: 

C 

C  THIS  SUBPROGRAM  COMPUTES  THE  GRID  VOLUME  LIQUID  WATER 
C  CONTENT,  QL,  AND  THE  FRACTIONAL  CLOUDINESS,  RC.  ALSO, 

C  THE  EDDY  TRANSFER  COEFFICIENTS  ARE  COMPUTED  —  EITHER 

C  USING  YAMADA'S  OR  THE  HASSID  AND  GALPERIN  (1983)  METHOD. 

C  AT  LAND  POINTS,  TIME-DEPENDENT  BUDGET  EQUATIONS  FOR  SURFACE 
C  TEMPERATURE  AND  MOISTURE  ARE  SOLVED  USING  THE  FORCE- 
C  RESTORE  METHOD.  SOMBER  I A  AND  DEARDORFF  (1977),  HASSID 
C  AND  GALPERIN  (1983),  AND  DEARDORFF  (1978)  ARE  PRIMARY 

C  SOURCES. 

C 

C  ORIGINAL  PROGRAM^R,  DATE:  STEPHEN  BURK  (NOARL),  1987 

C 

C  CURRENT  PROGRAhMER:  ROLF  LANGLAND  (NOARL) 


C  HOWARD  LEW  IT  (CSC) 

C  LIANA  ZAr«RESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  HOSKEEP 
C 

C . MAINTENANCE  SECTION. 

C 


C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  RECORD  OF  CHANKSS: 

C 

C . END  PROLOGUE 
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SUBROUTINE  INITIAL 

. START  PROLOGUE . 

SUBPROGRAM  NAME:  INITIAL 

DESCRIPTION: 

THIS  IS  THE  INITIALIZATION  SLBROUTINE  IN  WHICH  THE  NOQAPS 
OR  NORAPS  PROFILES  ARE  READ  IN  AND  INTERPOLATED  TO  THE  NABL 
GRID.  MANY  OF  THE  MODEL  ARRAYS  FOR  TURBULENCE  VARIABLES 
ARE  INITIALIZED  AS  ZERO  HERE.  ALSO,  INITIAL  VALUES  FOR 
RADIATIVE  FLUX  CALCULATIONS  ABOVE  THE  TOP  OF  THE 
COMPUTATIONAL  GRID  ARE  SET  HERE. 

ORIGINAL  PROGRAM-ER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAft^R:  ROLF  LANO-AND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  ZAMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  INITIAL 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 


PARAMETER  (NN=40,NR=21 ,NL=1 1 ,NT=37,NP=20) 

CHARACTER  CMOCH'S 

DIMENSION  RHOAU(NN),RHOVU(NN) 

C0^t10N/BLK1/ALF,VISK.A1,A2,B1,B2,C,G,BETA,VK,SQ,Z(NN), 

1  BT,R,QAW,GAMH,HEATL,CPO,RD,QRD,EFLUX(NN),SFLUX(NN) 

COM-ION/BLK2/PLOC(  10)  ,TTI  (NR),QQI  (NR)  ,UUI  (NR)  ,VVI  (NR)  ,PP  I  (NR)  , 

1  KBLKST,KTOT,CMOO 

OOm)N/BLANK/ABT  (  NN  ) ,  ABTH  (NN) ,  Cl  (  NN  ) ,  PS ,  ENUL  (  NN  )  ,  LL , 

2  C2(NN)  ,C3(NN)  ,CLAM2(NN)  ,UN(WI)  ,VN(NN)  ,UF(NN)  ,VF(NN)  , 

3  DQWDZ(NN) ,DT,DTHLDZ(NN) ,DUDZ(NN),DUMV(NN) ,DVDZ(NN) ,EKH(NN) , 

4  EKM(NN) , IM1 , IM2, IM3, IT, (TOT,P(NN) ,X1U(NN) ,X2U(NN) ,DTHVDZ(NN) , 

5  PRATIO(NN),QL(NN) ,QRN(NN) ,QSL(NN) ,QS(NN) ,X2UL(NN) , 

6  QW(NN),QWSB(NN),QWSTAR,QWZ,RC(NN),RLAT,RLON, IBLKST(NP), 

7  SPECHUM ( NN ) , PH  I M , OOSALP , S I NALP , SST , WK ( NN , NL ) , 

8  TH(NN) ,THL(NN) ,THLQWB(NN) ,TH_S8(NN) ,HDUCT,GAM2, 

9  THLSTAR,THLWB(NN) ,THLZ,THVWB(NN) ,TKE(NN) ,TS IG(NR) , 

A  TKESR(NN) ,TL(NN) ,TV(NN) ,U(NN) ,UG(NN) ,USTAR, V(NN) ,VG(NN) , 

B  WQWB(NN) ,XL(NN) ,THTEND(NN) ,QTEND(NN) ,VT(NN) ,TDEW(NN) , 
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C  Y1 ,Z1 ,ZL,ZN(NN) ,ZS(NN) ,DECLN,ZTIME,UTEND(NN) ,VTEND(NN) .PSTEND, 
D  QRL(NN),QRS(NN), TIMER, IDELAY,XMZ,XBC,ZS1G(NR) ,XM(NN) , 

E  m (NL) ,QQ1 (NL) ,FML(NN) ,FPL(NN) ,TSTEND ,GW,GWWG,TOPO, 

F  TT2(NL),QQ2(NL),FMS(NN),FPS(NN), 

G  SIQDTUNL),SIQDT2(NL) 

C 

C  THIS  S/R  READS  THE  INITIAL  CONDITIONS  COMING  FROM  NORAPS  OR  NOGAPS. 
C  THESE  INITIAL  CONDITIONS  AT  GRIDPOINT  JL  ARE: 

C 

C  PLOC(1)  =  IPT 

C  PLOC(2)  =  JPT 

C  PLOC(3)  =  LATITUDE 

C  PLOC(4)  =  LONGITUDE 

C  PLOC(5)  =  I  SEA  (CURRENTLY  NOT  USED) 

C  PLOC(6)  =  TSEA 

C  PLOC(7)  =  PS2 

C  PLOC(8)  =  GRDWET 

C  PLOC(9)  =  TOPO 

C  PLOC(IO)  =  TSINIT:  NORAPS  DEEP  SOIL  TEMP  (CURRENTLY  NOT  OUTPUT 
C  BY  OPERATIONAL  NORAPS) 

C  TTI(K)  =  TEMP  AT  SIGMA  LEVEL 

C  (HERE  K=1  EQUALS  K=KK  IN  NORAPS,  ETC.) 

C  QQI(K)  =  SPECIFIC  HUMIDITY 

C  UUI(K)  =  U-COMPONENT  WIND 

C  WI(K)  =  V-COMPONENT  WIND 

C  PPI(K)  =  INITIAL  PRESSURE  ON  SIGMA  LEVELS 

C 
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SUBROUTINE  NABL42( JTAU.NTAU) 
. START  PROLOGUE 


SUBPROGAM  NAME :  NABL42 

DESCRIPTION:  WRITE  NABL  FORECAST  TO  FILE  (FNOC  FORMAT)  FOR  SKEW-T 

ORIGINAL  PROGRAI-MER:  ROLF  LANGLAND  (1991) 

CURRENT  PROGRAMtR:  ROLF  LANGLAND  (NOARL) 

LIANA  ZAMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  NABL42( JTAU.NTAU) 

SET  NABL  NAMELIST  VARIABLE  "FNSKEW"  =  .TRUE. 

INPUT  PARAMETERS: 

JTAU  -  CORRESPONDS  TO  ITAU(M)  NABL  FORECAST  TAU 
NTAU  -  YYMIJOHH  OF  NABL  FORECAST  TAU  (CH*8) 

OUTPUT  PARAf«TERS: 

RLAT  -  LATITUDE  (0-90),  INS  (N  OR  S) 

BLON  -  LONGITUDE  (0-180),  lEW  (E  OR  W) 

CMOO  -  NORAPS  OR  NOGAPS  IDENTIFIER 

P  -  PRESSURE  (MB) 

T99  -  TEMPERATURE  (C) 

TD99  -  DEW  POINT  TEMPERATURE  (C) 

IV90  -  WIND  DIRECTION  IND  DEG  (0-360) 

1V99  -  WIND  SPEED  (KNOTS) 

ZFT  -  NABL  GRID  POINT  HEIGHT,  INC  TERRAIN  (FT) 

TOPO  -  TERRAIN  HEIGHT  (FT) 

OUTPUT  FILES:  TAPE42 


MAINTENENCE  SECTION 


RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  NABL46( JTAU.NTAU) 

. START  PROLOGUE . 

SUBPROGRAM  NAME:  NABL46 

DESCRIPTION:  WRITE  REFRACTIVITY  DATA  TO  FILE  IN  FNOC  FORMAT 


ORIGINAL  PROGRATTER:  ROLF  LANGLAND 

CWIRENT  PROGRATt^R:  ROLF  LANGLAND 

LIANA  2AMBRESKY 

USAGE  (CALLING  SEQUENCE):  CALL  NABL46( JTAU,NTAU) 
SET  NABL  NAMELIST  VARIABLE  •'RNRERC"  =  .TRUE. 


INPUT  PARAMETERS: 

JTAU  -  CORRESPONDS  TO  ITAU(M)  NABL  FORECAST  TAU 
NTAU  -  YYM130HH  OF  NABL  FORECAST  TAU  (CH*8) 


OUTPUT  PARAMETERS: 

RLAT  -  LATITUDE  (0-90),  INS  (N  OR  S) 

BLON  -  LONGITUDE  (0-180),  lEW  (E  OR  W) 

CMOD  -  NORAPS  OR  NOGAPS  IDENTIFIER 
P  -  PRESSURE  (MB) 

T99  -  TEMPERATURE  (C) 

TD99  -  DEW  POINT  TEMPERATURE  (C) 

XMZ  -  SURFACE  MODIFIED  REFRACTIVITY 

TOPO  -  TERRAIN  HEIGHT  (FT) 

XM  -  MODIFIED  REFRACTIVITY 

ZMDEL  -  LAPSE  RATE  OF  REFRACTIVITY 
RFRAC  -  REFRACTIVITY  CATEGORY 

ZFT  -  NABL  GRID  POINT  HEIGHT,  INC  TERRAIN  (FT) 


OUTPUT  FILES:  TAPE46 


MAINTENENCE  SECTION 


RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  NABL75( ITIME.NDTG, IMOO) 

. START  PROLOGUE . 

SUBPROGRAM  NAME;  NABL75 

DESCRIPTION:  WRITE  NABL  FORECAST  TO  FILE  FOR  SKEW-T  PLOT 

ORIGINAL  PROGRATtlER:  ROLF  LANGLAND  (1991) 

OURRENT  PROGRAM1ER:  ROLF  LANGLAND  (NOARL) 

LIANA  ZAMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  NABL75( ITIME.NDTG, IMOD) 

SET  NABL  NA^ELIST  VARIABLE  "NBSKEW"  =  .TRUE. 

INPUT  PARAMETERS: 

ITIME  -  NABL  FORECAST  TAU  (INTEGER) 

NDTG  -  YYMIXDHH  OF  NABL  FCST  TAU  (18) 

IMOD  -  1=NABL/NORAPS,  2=NABL/NOGAPS 

OUTPUT  PARAMETERS: 

KBLKST  -  BLOOKSTATION  IDENTIFIER  (5-DIGIT)  FOR  FORECAST  POINT 
PLOC(1)-  NORAPS/NOGAPS  GRID  POINT  IN  EAST-WEST  DIRECTION 
PLOC(2)-  NORAPS/NOGAPS  GRID  POINT  IN  NORTH-SOUTH  DIRECTION 
PLOC(5)-  NORAPS/NOGAPS  LAND-SEA- ICE  CODE 
PLOC(9)-  TERRAIN  HEIGHT  (M) 

RLAT  -  LATITUDE  (-90  S-HEM  TO  +90  N-HEM)  OF  FORECAST  POINT 
DLON  -  LONGITUDE  (-180  E-HEM  TO  +180  W-HEM)  OF  FORECAST  POINT 
OmQ  -  GROUND  MOISTURE  PARA^€TER 

SST  -  SEA  SURFACE  TEMPERATURE,  OR  GROUND  SKIN  TEMPERATURE 
P  -  PRESSURE  (MB) 

T99  -  TEMPERATURE  (C) 

TD99  -  DEW  POINT  TEMPERATURE  (C) 

UV9D  -  WIND  DIRECTION  IN  DEG  (0-360) 

UV99  -  WIND  SPEED  (M/S) 

IH  -  NABL  GRID  POINT  HEIGHT,  INC  TERRAIN  (M) 

OUTPUT  FILES:  TAPE75 


MAINTENENCE  SECTION 


RECORD  OF  CHANGES 


END  PROLOGUE 
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SUBROUTINE  NABL80(  ITIME.NDTG.X.NK) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  NABL80 
C 

C  DESCRIPTION:  WRITE  NORAPS  FORECAST  TO  FILE  FOR  SKEW-T  PLOT 
C 

C  ORIGINAL  PROGRAhTER:  ROLF  LANGLAND  (1991) 

C 

C  CURRENT  PROQRArt^R:  ROLF  LANGLAND  (NOARL) 


C  LIANA  ZAMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  NABL80( ITIME.NDTG, JL.NK) 

C  SET  NABL  NAMELIST  VARIABLE  "NRSKEW"  =  .TRUE. 

C 

C  INPUT  PARAMETERS: 

C  I  TIME  -  NABL  FORECAST  TAU  (INTEGER) 

C  NDTG  -  YVmXXH  OF  NABL  FORECAST  TAU  (18) 

C  JL  -  COUNTER  FOR  SEQUENCE  OF  NABL  FORECAST  POINTS 

C  NK  -  COUNTER  FOR  SEQUENCE  OF  NORAPS  SKEW-T  SAVE  TIMES 

C 

C  OUTPUT  PARAMETERS: 

C  KBLKST  -  BLOCKSTATION  IDENTIFIER  (5-DIGIT)  FOR  FORECAST  POINT 

C  RLAT  -  LATITUDE  (-90  S-HEM  TO  +90  N-HEM)  OF  FORECAST  POINT 

C  DLON  -  LONGITUDE  (-180  E-HEM  TO  +180  W-HEM)  OF  FORECAST  POINT 

C  IMOO  -  (  =  4)  IDENTIFIES  NORAPS  IN  SKEW-T  SOFTWARE 

C  SLOC3(1,.)  -  NORAPS  GRID  POINT  IN  X-DIRECTION 

C  SLOC3(2,.)  -  NORAPS  GRID  POINT  IN  Y-DIRECTION 

C  PLOC(8)-  GROUND  MOISTURE  PARAMETER 

C  PLOC(9)-  TERRAIN  HEIGHT  (M) 

C  PLOC(5)-  NORAPS  LAND-SEA- ICE  CODE 

C  SSTNR  -  SEA  SURFACE  TEMPERATURE  OR  GROUND  SKIN  TEMPERATURE 

C  PPS  -  PRESSURE  (MB) 

C  T9 9  -  TEMPERATURE  (C) 

C  TD99  -  DEW  POINT  TEMPERATURE  (C) 

C  UV9D  -  WIND  DIRECTION  IN  DEG  (0-360) 

C  UV99  -  WIND  SPEED  (M/S) 

C  NRZ  -  HEIGHT  OF  NORAPS  SIGMA  LEVELS  (M) 

C 

C  OUTPUT  FILES:  TAPE80 
C 

C . MAINTENENCE  SECTION . 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 


C-I7 


ooooooooooooooooooooooooooooooooo 


NABLS  SDD/IDD  Original 


1 5  Noveaabe  r  19  91 


SUBROUTINE  OUT( IL, JL,CLAT,CLON,DSKIP) 

. START  PROLOGUE . 

SUBPROGRAM  NAME :  OUT 

DESCRIPTION: 

TVIIS  IS  THE  MAIN  OUTPUT  SUBROUTINE  FOR  NABL.  LOGICAL 
VARIABLE,  DSKIP,  IS  USED  TO  CONTROL  THE  VOLUME  OF 
OUTPUT  ONE  GETS. 

ORIGINAL  PROGRAI-MER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAMMER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  2AMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  OUT( I L , JL ,CLAT,CLON,DSK I P) 

INPUT  PARAMETERS: 

1 L , JL  =  NOGAPS/NORAPS  I , J  VALUES  FOR  SELECTED  NABL  PO 1  NT 
CLAT,CLON  =  NOGAPS/NORAPS  LATITUDE  &  LONGITUDE  OF  NABL  POINT 
DSKIP  =  LOGICAL  VARIABLE  USED  TO  DETERMINE  WHETHER  TO  SKIP 
PRINTING  OUT  A  BLOCK  OF  INFORMATION 

. MAINTENANCE  SECTION . 

PRINIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  TKDAG(NDTG,MT) 

C 

C  SUBPROGRAM  NAME:  TMDAG 
C 

0 . START  PROLOGUE . 

C 

0  DESCRIPTION:  WRITE  SELECTED  NABL  FORECAST  VARIABLES  TO  OUTPUT 
C  FILE  FOR  TIME  SERIES  PLOTTING 

0 

0  ORIGINAL  PROGRATTER:  ROLF  LANGLAND 
C 

C  CURRENT  PROGRAMER:  ROLF  LANGLAND 

0  LIANA  ZAMBRESKY 

0 

C  USAGE  (CALLING  SEQUENCE):  CALL  TMDAG(NDTG,MT) 

C  SET  NABL  NAMELIST  VARIABLE  "TMSER’*  =  .TRUE. 

C 

0  INPUT  PARAMETERS: 

0  NDTG  -  YYMMDDHH  OF  NABL  'FORECAST  TAU  (l8) 

C  MT  -  INTEGER  FLAG  FOR  TAU-0  OR  NOT  TAU-0 

C 

0  OUTPUT  PARAMETERS: 

C  KBLKST  -  BLOCKSTATION  IDENTIFIER  (5-DIGlT)  OF  FORECAST  POINT 

C  THL  -  LIQUID  WATER  POTENTIAL  TEMPERATURE  (K) 

0  QW  -  TOTAL  WATER  CONTENT  (KG/KG) 

0  UV99  -  WIND  SPEED  (M/S) 

C  UV9D  -  WIND  DIRECTION  (0-360  DEG) 

0  TKE  -  TURBULENT  KINETIC  ENERGY 

C  THVWB  -  TURBULENT  SENSIBLE  HEAT  FLUX  (WATTS  PER  SQUARE  METER) 

C  WQWB  -  TURBULENT  LATENT  HEAT  FLUX  (WATTS  PER  SQUARE  METER) 

C 

C  OUTPUT  FILES:  TAPE70  (ALIAS  PDATA) 

C 

0 . MAINTENENCE  SECTION . 

0 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  TQTEND(IHR) 

. START  PROLOGUE . 

SUBPROGRAM  NAME :  TQTEND 

DESCRIPTION: 

THIS  SUBPROGRAM  COMPUTES  THE  HOURLY  AVERAGE  ADIABATIC 
TENDENCIES  FOR  TEMPERATURE  AND  MOISTURE  BASED  ON  NOGAPS/ 
NORAPS  OUTPUT.  CODE  WITHIN  NOGAPS/NORAPS  SAVES  THE 
TEMPERATURES  AND  MOISTURES  THAT  RESULT  AT  HOURLY  INTERVALS 
WHEN  ONLY  ADIABATIC  (ADVECTIVE)  PROCESSES  ARE  CONSIDERED. 
THUS,  THESE  TENDENCIES,  ONCE  INTERPOLATED  TO  THE  NABL 
GRID,  BECOME  THE  ADIABATIC,  ADVECTIVE  FORCING  THAT  GIVES 
NABL  ITS  3-D  CHARACTER  I.E.,  THEY  PERMIT  NABL  TO 
RESPOND  TO  SYNOPTIC  CHANGES  AND  SUBSIDENCE  IN  A  MANNER 
THAT  A  STRICTLY  1-D  MODEL  WOULD  NOT. 

ORIGINAL  PROGRAMMER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAMMER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  2AMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE):  CALL  TQTEND(IHR) 

INPUT  PARAMETERS: 

I  HR  =  FORECAST  HOUR  AT  WHICH  THIS  SUBROUTINE  IS  CALLED 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 

. END  PROLOGUE . 
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SUBROUTINE  UTENDS(IHR) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME;  UTENDS 
C 

C  DESCRIPTION: 

C 

C  THIS  SUBPROGRAM  COMPUTES  THE  WIND  TENDENCIES  BASED  ON 
C  THE  NOGAPS/NORAPS  WIND  COMPONENTS  AT  ONE  HOUR  INTERVALS. 

C  THESE  TENDENCIES  ARE  USED  TO  NUDGE  THE  NABL  WIND  FORECAST 
C  TOWARDS  THE  NOGAPS/NORAPS  FORECAST  SO  THAT  IN  THE  ABSENCE 
C  OF  TURBULENT  MIXING  THE  NABL  FORECAST  WINDS  ARE  VERY 
C  SIMILAR  TO  THE  NOGAPS/NORAPS  WINDS.  THE  TENDENCIES  ARE 
C  ALSO  INTERPOLATED  TO  THE  NABL  GRID  HERE. 

C 

C  ORIGINAL  PROGRAMMER ,  DATE:  STEPHEN  BURK  (NOARL),  1987 

C 

C  CURRENT  PROGRAMMER:  ROLF  LANGLAND  (NOARL) 

C  HOWARD  LEW  IT  (CSC) 

C  LIANA  ZAMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  UTENDS (1 HR) 

C 

C  INPUT  PARAMETERS: 

C 

C  IHR  =  FORECAST  HOUR  AT  WHICH  THIS  SUBROUTINE  IS  CALLED 
C 

C . MAINTENANCE  SECTION . 


C 

C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE 

C 
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SUBROUTINE  UVCOMP(PHI ,PHIN,EKM,ZS,ZN,WK,DT,USTAR,ANG,Z) 
. START  PROLOGUE . 


SUBPROGRAM  NAME:  UVCOMP 

DESCRIPTION: 

THIS  SUBPROGRAM  COMPUTES  THE  MOMENTUM  FLUX  DIVERGENCE 
FOR  A  GIVEN  WIND  OOMPONENT.  A  TRIDI AGONAL  IMPLICIT 
EC3UATION  IS  SOLVED  FOR  THIS  EDDY  MIXING  PROCESS.  THE 
SOLUTION  IS  NUDGED  TOWARDS  THE  NOGAPS/NORAPS  WINDS  — 

THIS  NUDGING  BEING  DONE  TOWARDS  THE  END  OF  THE  MAIN 
PROGRAM.  IN  REGIONS  WHERE  THE  MIXING  IS  WEAK,  THE 
NUDGING  WILL  CAUSE  THE  NABL  WIND  PROFILE  TO  EVOLVE 
SIMILAR  TO  NOGAPS/NORAPS. 

ORIGINAL  PROGRAWER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

CURRENT  PROGRAI-MER:  ROLF  LANGLAND  (NOARL) 

HOWARD  LEW  IT  (CSC) 

LIANA  ZAMBRESKY  (FNOC) 

USAGE  (CALLING  SEQUENCE): 

CALL  UVCOMP(PHI ,PHIN,EKM,ZS,ZN,WK,DT,USTAR,ANG,Z) 

INPUT  PARAMETERS: 

PHI  =  WIND  OOMPONENT  (U  OR  V)  AT  TIME  LEVEL  N 
EKM  =  EDDY  COEFFICIENT  FOR  MOMENTUM 

ZS  =  GRID  SPACING  DIFFERENCE  BETWEEN  CURRENT  LEVEL,  1,  AND  1-1. 

ZN  =  GRID  SPACING  DIFFERENCE  BETWEEN  CURRENT  LEVEL,  I,  AND  1+1. 

WK  =  GENERAL  PURPOSE  WORK  ARRAY 

DT  =  TIME  STEP  (SEC) 

USTAR  =  FRICTION  VELOCITY  (  M/S  ) 

ANG  =  COSINE  OR  SINE  OF  SURFACE  WIND  "^0  »0S1TIVE  X-AXIS 

Z  =  GRID  POINT  HEIGHT  ARRAY 

OUTPUT  PARAMETERS: 

PH  IN  =  VALUE  OF  WIND  COMPONENT  (U  OR  V)  AT  NEW  TIME  LEVEL,  N+1. 

. MAINTENANCE  SECTION . 

PRINCIPLE  VARIABLES  AND  ARRAYS: 

RECORD  OF  CHANGES: 


END  PROLOGUE 
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SUBROUTINE  UVSD  (U,V,SPD,DIR,N) 

C 

C . START  PROLOGUE 

C 

C  SUBPROGRAM  NAME:  UVSD 
C 

C  DESCRIPTION: 

C 


C  THIS  SUBPROGRAM  TAKES  U  AND  V  WIND  COMPONENTS  AND  CONVERTS 
C  THEM  TO  SPEED  (M/S)  AND  DIRECTION  (DEG). 

C 

C  ORIGINAL  PROGRAI-MER,  DATE:  ROLF  LANGLAND  (NOARL),  1990 

C 

C  CURRENT  PROGRATtlER:  ROLF  LANGLAND  (NOARL) 

C  HOWARD  LEW IT  (CSC) 

C  LIANA  ZAMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  UVSD(U,V,SPD,D1R) 

C 

C . MAINTENANCE  SECTION . 

C 

C  INPUT  PARAMETERS: 

C 

C  U  =  EAST-WEST  WIND  COMPONENT  (  M/S  ) 

C  V  =  NORTH-SOUTH  WIND  COMPONENT  (  M/S  ) 

C 

C  OUTPUT  PARAMETERS: 


C 

C  SPD  =  TOTAL  WIND  SPEED  (  M/S  ) 

C  DIR  =  WIND  DIRECTION  FROM  NORTH  (  DEG  ) 
C 

C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE. . . 
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SUBROUTINE  VINTRP(X,Z,XSD, J1 ,K1 , lEXTRP) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  VINTRP 

C 

C  DESCRIPTION: 

C 

C  THIS  SUBPROGRAM  INTERPOLATES  NOGAPS  OR  NORAPS  VALUES  TO 
C  1>C  NABL  GRID.  THE  INTERPOLATION  IS  DONE  LINEARLY  IN 

C  HEIGHT.  IF  IEXTRP=1,  NABL  POINTS  BELOW  THE  LOWEST  NOGAPS/ 

C  NORAPS  POINT  ARE  SET  EQUAL  TO  THE  NOGAPS/NORAPS  VALUE  AT 

C  THEIR  LOWEST  GRID  POINT  ABOVE  THE  SURFACE.  IF  lEXTRP  IS 

C  NOT  EQUAL  TO  ONE,  THEN  NABL  POINTS  BELOW  THE  LOWEST  NOGAPS/ 

C  NORAPS  POINT  ARE  FOUND  BY  EXTRAPOLATING  THE  NEAR  SURFACE 
C  NOGAPS/NORAPS  GRADIENT  DOWNWARD. 

C 

C  ORIGINAL  PROGRATtlER,  DATE:  STEPHEN  BURK  (NOARL),  1987 

C 

C  CURRENT  PROGRATTIER:  ROLF  LANGLAND  (NOARL) 

C  HOWARD  LEW IT  (CSC) 

C  LIANA  ZAMBRESKY  (FNOC) 

C 

C  USAGE  (CALLING  SEQUENCE):  CALL  V1NTRF(X,Z,XSD, J1 ,K1 , lEXTRP) 

C 

C  INPUT  PARAMETERS; 

C 

C  Z  =  NABL  GRID  POINT  HEIGHTS 

C  XSD  =  VARIABLE  VALUE  AT  NOGAPS/NORAPS  LEVEL 

C  J1,K1  =  VALUES  USED  TO  SET  SECOND  INDEX  OF  ARRAY  XSD 

C  lEXTRP  =  DETERMINES  NATURE  OF  INTERPOLATION/  EXTRAPOLATION 
C  FOR  NABL  POINTS  LOWER  THAN  LOWEST  NOGAPS/NORAPS  PT. 

C 

C  OUTPUT  PARAMETERS : 

C 

C  X  =  INTERPOLATED  VALUE  ON  NABL  GRID  OF  THE  INPUT  VARIABLE,  XSD. 
C 

C . MAINTENANCE  SECTION . 

C 

C  PRINCIPLE  VARIABLES  AND  ARRAYS: 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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PROGRAM  NGNPNT( INPUT, OUTPUT, TER, TAPE1= I NPUT.TAPE 2 =TER, 

♦  PRO, TAPE7=PRO, TAPES) 

C 

C . START  PROLOGUE . 

C 

C  PROGRAM  NAME:  PROGRAM  NGNPNT 
C 

C  DESCRIPTION: 

C  INTERACTIVE  ENTRY  OF  NOGAPS  NABL  POINTS  FOR  TRANSFER  TO 
C  DGFILE  ON  SAM 
C 

C  ORIGINAL  PROGRAI^R,  DATE:  HOWARD  LEWIT,  20  DEC  1990 
C 

C  CURRENT  PROGRAWER: 

C 

C  CONTRACT  NUMBER  AND  TITLE:  NONE 
C 

C  REFERENCES:  NONE 
C 

C  COMPUTER/OPERATING  SYSTEM:  FNOC  NOS  BE  PEPS 
C 

C  LIBRARIES  OF  RESIDENCE:0PSPL1 
C 

C  CLASS F  1  CAT !  ON :  UNCLAS 
C 

C  USAGE  (XL): 

C  EXECUTED  IN  INTERACTIVE  PROCEDURE  NQNABL: 

C  APLIB(MT1731,*NGNPNT) 

C  NGNPNT. 

C  (SEE  PROCEDURE  NGPO I  NT , OPSF ! L  FOR  TREATMENT  OF  OUTPUT) 

C  INPUT  FILES:  TAPE1=  INTERACTIVE  USER  INPUT 
C 

C  OUTPUT  FILES: 

C  TAPE2  =  DIRECTIVES  TO  INTERACTIVE  USER 
C  TAPES  =  DATA  FOR  TRANSFER  TO  SAM 
C 

C  ERROR  CONDITIONS: 

C  NUMEROUS  MESSAGES  TO  USER  FOR  INCORRECT  ENTRIES  AND  FORMATS 
C 

C  ADDITIONAL  OOTtlENTS: 

C 

C . MAINTENANCE  SECTION . 

C 

C  BRIEF  DESCRIPTION  OF  PROGRAM  MODULES: 

C  SUBROUTINE  NQOLD  (  &  ENTRY  NQWRT)  READ  AND  WRITE  THE 
C  EXISTING  POINT  FROM  FILE  NQNREC 
C  SUBROUTINE  NGRPT  -  USER  INPUT  NGPO I NT  LOCATION 
C  SUBROUTINE  NGFORM  -  CHANGES  THE  POINTS  TO  SAM  INPUT  FORMAT  FOR 
C  TAPE  S 

C 

C  DATA  FILES:  FILE  NGNREC  CONTAINS  POINTS  ON  FRONT-END  FILE 
C  FOR  REFERENCE  AND  BACKUP 

C 

C  TEMPORARY  FILES:  NONE 
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C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS:  SEE  COMDECK  D I  SCRIPT  IONS 
C 

COMMON  BLOCKS:  SEE  COMDECKS 


METHOD: 

DISPLAYS  EXISTING  POINTS  TO  USER,  PROMPTS  CHANGES  TO  INDIVIDUAL 
POINTS.  USER  ENTERS  POINTS  IN  STANDARD  FORMAT  ,  PROGRAM  CHANGES 
THEM  TO  NOGAPS  REQUIRED  FORMAT  AND  WRITES  TO  TAPES 
C  PROCEDURE  NGNMOD  TRANSFERS  TAPES  TO  SAM  AND  WRITES  TO 
C  DGFILE. 

C 

C  LANGUAGE:  FTN5 
C 

C  RECORD  OF  CHANGES: 

C  «CHANGE  NOTICE>>  NGNPNT*01  (9  JAN  1991)  --  LEWIT 
C  MODIFY  COM-ENTS 

C . END  PROLOGUE . 
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SUBROUTINE  NGEXTR(L I NE ,  I  S ,  IN,  IT,FF,NUM,RR,NX,=*') 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NGEXTR 
C 

C  DESCRIPTION: 

C  EXTRACT  CHARACTERS,  INTEGERS  OR  REAL  NUMBERS  FROM  CHARACTER 
C  STRING 
C 

C  ORIGINAL  PROGRAM-ER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGEXTR(LINE, IS, IN, IT,FF,NUM,RR,NX,*1) 

C  ALTERNATE  RETURN  *1  FOR  BLANK  INPUT 
C 

C  INPUT  PARAMETERS: 

C  LINE  -  80  CHARACTER  VARIABLE 
C  IS  -  STARTING  POINT  IN  LINE 
C  IN  -  MAXIMUM  NUMBER  OF  CHARACTERS  TO  SEARCH 
C  IT  -  OUTPUT  TYPE 
C  0  -  CHARACTER 

C  1  -  INTEGER 

C  2  -  REAL  number 

C 

C  OUTPUT  PARAMETERS: 

C  FF  -  CHARACTER  STRING 
C  NUM  -  INTEGER 
C  RR  -  REAL  NUMBER 
C 

C  INPUT  FILES:  NONE 

C 

C  OUTPUT  FILES:  NONE 
C 

C  OOMiON  BLOCKS:  NONE 
C 

C  ERROR  CONDITIONS:  NONE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  SEE  I NPUT/OLTTPUT  PARAMETERS 
C  FIRST  -  FLAG  FOR  FIRST  UNIT 
C  NCHR  -  NUMBER  OF  UNIT  COUNT 
C  DEC  -  LOCATION  IN  NCHR  OF  DECIMAL  TYPE  2 
C  KNA/KNB  -  ASCI  I  NUMBER  RANGE 
C  KCA/KCB  -  ASCI  I  LETTER  RANGE 
C  NSGN  -  (')  SIGN  FLAG 
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C 

C  METHOD: 

C  IDENTIFY  EACH  CHARACTER  BY  TYPE  AND  ACCUMULATE  THEM 
C  ACCORDING  TO  DECODE  TYPE. 

C  PROVIDE  (-)  IF  NOT  CHARACTER  AND  (.)  IF  TYPE  2 
C 

C  RECORD  OF  CHANGES: 

C  <<CHANGE  NOTICE>>  NGEXTR*01  (20  MAR  91)  —  LEWIT 
C  FIX  DECIMAL  PLACE  FOR  REAL  NUMBERS 
C 

C . END  PROLOGUE . 
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SUBROUTINE  NGFORM 
C 

C . START  PROLOGUE . 

0 

C  SUBPROGRAM  NAME:  SUBROUTINE  NGFORM 
0 

C  DESCRIPTION: 

0  CONVERT  TO  NOGAPS  INPUT  FORM  AND  WRITE  OUT  IN  NAI^LIST  FORMAT 
C 

C  ORIGINAL  PROGRAM«R,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUf^ER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGFORM 
C 

C  1 NPUT  PARAMETERS : 

C  ARRAYS  OF  LAT,  LONG  AND  LAND/OCEAN 
C 

C  OUTPUT  PARAMETERS: 

C  rWNE 

c 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  TAPE  8 
C 

C  O0^t1ON  BLOCKS:  $NGPNT 
C 

C  ERROR  CONDITIONS:  BAD  NAMELIST  WRITE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  SEE  INPUT/OUTPUT  PARAMETERS 
C 

C  METHOD: 

C  CONVERT  LAT,  LONG  TO  NOGAPS  INPUT  FORMAT 
C  LAT  NORTH  =  +,  SOUTH  =  - 
C  LONG  E  =  +,  WEST  =  - 

C  WRITE  NAMELIST  TO  TAPE  8  FOR  TRANSFER  TO  SAM 
C 

C  RECORD  OF  CHANGES: 

C 

C  <<CHANGE  N0TICE>>  NGFORM*01  (20  MAR  1991)  —  LEWIT 
C  INCREASE  LAT/LON  UP  TO  3  DECIMAL  PLACES 

C . END  PROLOGUE . 
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SUBROUTINE  NGLL(AL ,NS, AN ,EW,0) 

C 

C. . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NGLL 
0 

C  DESCRIPTION: 

C  PROVIDE  FOR  INPUT  OF  LAT/LONG  POSITIONS 
C 

C  ORIGINAL  PROGRAMMER,  DATE:  HOWARD  LEW IT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGLL(AL,NS,AN,EW,0) 

C 

C  INPUT  PARAMETERS: 

C  NONE 
C 

C  OUTPUT  PARAMETERS: 

C  AL.NS  =  LAT 
C  AN,EW  =  LONG 
C  0  =  LAND/OCEAN 
C 

C  INPUT  F1LES:TAPE1 
C 

C  OUTPUT  FILES:  TAPE 2 
C 

C  COrnON  blocks  : 

c 

C  ERROR  CONDITIONS: 

C  NUhCROUS  MESSAGES  TO  USERS  IN  CASE  OF  INVALID  INPUT 


C  UP  TO  NG  ERRORS  ALLOWED  THEN  MESSAGE  TO  GET  HELP  AND  STOP 
C 

C . MAINTENANCE  SECTION . 

C 


C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  LINE-  80  CHARACTER  USER  INPUT 
C  I  TRY  -  ERROR  COUNT 
C 

C  METHOD: 

C  PROMPT  USER  TO  ENTER  FULL  LAT  LONG.  USER  SUBROUTINE  NGEXTR  TO 
C  EXTRACT  THE  NUMBERS  AND  CHARACTERS  FROM  THE  ENTRY.  CHECK  THE 
C  ENTRY  FOR  VALID  VALUES.  VERIFY  WITH  USER. 

C  IF  NOT  OORRECT  TRY  AGAIN  NG  TIMES. 

C  USER  ENTER  "LAND"  OR  "OCEAN" 

C 

C  RECORD  OF  CHANGES: 

C  <'^CHANGE  NOTICE>>  NGLL*01  (20  MAR  91)  —  LEWIT 
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INCREASE  LAT/LON  UP  TO  3  DECIMAL  PLACES 
. END  PROLOGUE . 


C-31 


oooooooooooooooooooooooooooooooooooooooooooooooooooooo 


NABLS  SDD/IDD  Original 


15  November  1991 


SUBROUTINE  NQOLD 

. . START  PROLOGUE . 

SUBPROGRAM  NAME:  SUBROUTINE  NQOLD 
DESCRIPTION: 

READ  AND  WRITE  POINT  INFORMATION  IN  STANDARD  FORMAT  AS 
PERMANENT  FILE  NQNREC,  ID=OP,  SN=SHARDSK 

ORIGINAL  PROGRA^^ER,  DATE:  HOWARD  LEW  IT,  DECEMBER  1990 

CONTRACT  NUMBER  AND  TITLE: NONE 

CLASSIFICATION:  UNCLASS 

USAGE  (CALLING  SEQUENCE): 

CALL  NQOLD 

INPUT  PARAMETERS: 

ARRAYS  OF  LAT,  LONG  AND  LAND/OCEAN 

OUTPUT  PARAMETERS: 

NONE 

INPUT  FILES:  NONE 

OUTPUT  FILES:  ZRANDIO  RECORD  NGNREC 

OUTPUT  FILE:  PERMANENT  FILE  NGNREC  ID=OP  ,  SN=SHARDSK 

OOrnON  BLOCKS:  SNGPNT 

ERROR  CONDITIONS:  RECORD  NOT  AVAILABLE  CONTINUE  GOING 
. MAINTENANCE  SECTION . 


DATA  FILES:  NONE 

TEMPORARY  FILES:  NONE 

PRINCIPAL  VARIABLES  AND  ARRAYS: 

SEE  INPUT/OUTPUT  PARAMETERS 

r«THOO: 

READ  AND  WRITE  THE  DATA  AS  UNFORMATTED  RECORDS  TO 
TAPES 

CATALOG  NEW  CYCLE  OF  NQNREC 
PURGE  OLDEST  CYCLE  OF  NGNREC 

RECORD  OF  CHANGES: 

<<CHANGE  N0T1CE>>  NQOLD*01  (9  JAN  1991)  —  LEVJIT 
MODIFY  OOrtCNTS 

<<CHANGE  ,N0T:CE>>  NQOLD*02  (24  APRIL  1991)  —  LEW  IT 
LIMIT  NUMBER  OF  POINTS  READ  IN  TO  NPT 
. END  PROLOGUE . 
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SUBROUTINE  NGRPT 
C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NGRPT 
C 

C  DESCRirTlON: 

C  USER  CHANGE,  DELETE  OR  ADO  POINTS  TO  DGFILE 
C 

C  ORIGINAL  PROGRATtCR,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGRPT 
C 

C  INPUT  PARAMETERS:  NONE 
C 

C  OUTPUT  PARAMETERS:  LAT  AND  LONG  IN  ARRAYS  ALAT  AND  ALON 
C 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  TAPE2 
C 

C  COMMON  BLOCKS:  NGPNT 
C 

C  ERROR  CONDITIONS: 

C  LIMIT  OF  NPT  ENTRY  POINTS,  MESSAGE  AND  CONTINUE 
C  INCORRECT  FORMAT  OF  POINT 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEYPCRARY  =^!LES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  AL  -  USER  INPUT  LAT 

C  AN  -  USER  INPUT  LONG 

C  NS/EW  -  USER  INPUT  N  OR  S  AND  E  OR  W 
C  NMP  -  COUNT  OF  NUMBER  OF  INPUT  POINTS 
C  ALAT  -  LAT  ARRAY 

C  ALON  -  LONG  ARRAY 

C  LS  =  LAND/OCEAN  ARRAY 

C 

C  METHOD: 

C  GET  USER  INPUT  OF  LAT,  LONG  AND  LAND/SEA  VIA  SUBROUTINE  NGLL 
C  SAVE  MODIFIED  OR  ADDED  DATA 
C  TEST  FOR  ADDITIONAL  INPUT  POINTS  UP  TO  NPT 
C  DELETE  INDICATED  POINT,  CLOSE  SPACES  IN  ARRAY 
C 

C  RECORD  OF  CHANGES: 

C 
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C  <<CHANGE  NOTICE>>  NGRPT*01  (9  JAN  1991)  —  LEWIT 
C  MODIFY  COIt^ENTS,  CORRECT  INITIAL  DATA  LIST  DISPLAY 
C<<CHANGE  NOTICE>>  NQRPT*02  (20  MAR  1991)—  LEWIT 
INCREASE  LAT/LON  UP  TO  3  DECIMAL  PLACES 
. END  PROLOGUE . 


SUBROUTINE  NGSTOP(IP) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  NGSTOP 
C 

C  DESCRIPTION: 

C  TERMINATE  THE  PROGRAM 
C 

C  ORIGINAL  PROGRAMMER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGSTOP(IP) 

C 

C  INPUT  PARAMETERS: 

C  IP  -  FLAG  1=  NORMAL,  2  =  ABNORMAL 
C 

C  OUTPUT  PARAMETERS:  NONE 
C 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  NONE 
C 

C  OOMMON  BLOCKS:  NONE 
C 

C  ERROR  CONDITIONS:  NONE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  IP  -  FLAG 
C 

C  METHOD: 

C  STOP  PROGRAM  WITH  MESSAGE  NORMAL  OR  ABNORMAL  TERMINATION 


C  DISCONNECT  TAPE1(  INPUT)  AND  TAPE2(  OUTPUT) 
C  STOP 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NGYN(*,*) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NGYN 
C 

C  DESCRIPTION: 

C  PROVIDE  PROMPT  TO  USER  TO  ENTER  Y  OR  N 
C 

C  ORIGINAL  PROGRAhMER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NGYN(*1,*2) 

C  ALTERNATE  RETURNS: 

C  *1  Y  RESPONSE 

C  *2  N  RESPONSE 

C 

C  INPUT  PARAMETERS:  NONE 
C 

C  OUTPUT  PARAMETERS:  NONE 
C 

C  INPUT  FILES:  TAPE1 
C 

C  OUTPUT  FILES:  TAPE 2 
C 

C  COmON  BLOCKS:  NONE 
C 

C  ERROR  CONDITIONS: 

C  USER  MUST  ANSWER  Y  OR  N,  KEEP  TRYING 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  YN  -  USER  RESPONSE 
C 

C  METHOD: 

C  PROMPT  USER  TO  ENTER  Y  OR  N.  IF  INCORRECT  KEEP  TRYING 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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PROGRAM  NRPPNT( INPUT, OUTPUT, TER, TAPE1=INPUT,TAPE2=TER, 

♦  PRO, TAPE7=PRO, TAPES) 

C 

C . START  PROLOGUE . 

C 

C  PROGRAM  NAME:  PROGRAM  NRPPNT 
0 

C  DESCRIPTION: 

C  INTERACTIVE  ENTRY  OF  NORAPS  NABL  POINTS  FOR  TRANSFER  TO 
C  NRNFIL7  ON  SAM 
C 

C  ORIGINAL  PROGRAMMER,  DATE:  HOWARD  LEWIT,  20  DEC  1990 
C 

C  CURRENT  PROGRAMMER: 

C 

C  CONTRACT  NUMBER  AND  TITLE:  NONE 
C 

C  REFERENCES:  NONE 
C 

C  COMPUTER/OPERATING  SYSTEM:  FNOC  NOS  BE  PEPS 
C 

C  LIBRARIES  OF  RESIDENCE:OPSPL1 
C 

C  CLASSIFICATION:  UNCLAS 
C 

C  USAGE  (XL): 

C  EXECUTED  IN  INTERACTIVE  PRXEDURE  NRPNABL: 

C  APLIB(Mm731, ♦NRPPNT) 

C  NRPPNT. 

C  (SEE  PROCEDURE  NRPNABL ,OPSF I L  FOR  TREATMENT  OF  OUTPUT) 

C 

C  INPUT  FILES:  TAPE1=  INTERACTIVE  USER  INPUT 
C 

C  OUTPUT  FILES: 

C  TAPE2  =  DIRECTIVES  TO  INTERACTIVE  USER 
C  TAPES  =  DATA  FOR  TRANSFER  TO  SAM 
C 

C  ERROR  CONDITIONS: 

C  NUMEROUS  MESSAGES  TO  USER  FOR  INCORRECT  ENTRIES  AND  FORMATS 
C 

C  ADDITIONAL  COMMENTS: 

C  NONE 

C . MAINTENANCE  SECTION . 

C 

C  BRIEF  DESCRIPTION  OF  PROGRAM  MODULES: 

C  SUBROUTINE  NROLD  -  READ  AND  WRITE  TWE  EXISTING  POINTS  FROM  PF 
C  NRNREC(AREA) 

C  SUBROUTINE  NRRPT  -  USER  INPUT  NRPNABL  LOCATION 
C  SUBROUTINE  NRFORM  -  CHANGES  THE  POINTS  TO  SAM  INPUT  FORMAT  FOR 
C  TAPE  8 

C  SUBROUTINE  NRPROC  -  PREPARE  PROCEDURE  PRO  TO  PASS  AREA  NUMBER  TO 
C  PROCEDURE  NRNMOO. 

C 

C  DATA  FILES:  FILE  NRPNABL  CONTAINS  FRONT  END  POINTS  FOR 
C  REFERENCE 
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C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS:  SEE  COMDECK  D I  SCRIPT  IONS 
C 

C  OOl-MON  BLOCKS:  SEE  OOMDECKS 
C 

C  METHOD: 

C  DISPLAYS  EXISTING  POINTS  TO  USER,  PROMPTS  CHANGES  TO  INDIVIDUAL 
C  POINTS.  USER  ENTERS  POINTS  IN  STANDARD  FORMAT  ,  PROGRAM  OHANGES 
C  THEM  TO  NORAPS  REQUIRED  FORMAT  AND  WRITES  TO  TAPES 
C  PROCEDURE  NRNMOO  TRANSFERS  TAPES  TO  SAM  AND  WRITES  TO 
C  NRNF I L7 . 

C 

C  LANGUAGE :  FTN5 
C 

C  RECORD  OF  CHANGES: 

C . END  PROLOGUE . 
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SUBROUTINE  NREXTR(L I NE , I S , IN, IT,FF,NUM,RR,NX,*) 

. START  PROLOGUE . 

SUBPROGRAM  NAME:  SUBROUTINE  NREXTR 
DESCRIPTION: 

EXTRACT  CHARACTERS,  INTEGERS  OR  REAL  NUMBERS  FROM  CHARACTER 
STRING 

ORIGINAL  PROGRA^^ER,  DATE:  HOWARD  LEW  IT,  DECEMBER  1990 

CONTRACT  NUMBER  AND  TITLE: NONE 

CLASSIFICATION:  UNCLASS 

USAGE  (CALLING  SEQUENCE): 

CALL  NREXTR(L1NE, I S, I N , IT ,FF ,NUM,RR,NX,*1 ) 

ALTERNATE  RETURN  *1  FOR  BLANK  INPUT 

INPUT  PARAMETERS: 

LINE  -  80  CHARACTER  VARIABLE 
IS  -  STARTING  POINT  IN  LINE 
IN  -  MAXIMUM  NUMBER  OF  CHARACTERS  TO  SEARCH 
IT  -  OUTPUT  TYPE 
0  -  CHARACTER 

1  -  INTEGER 

2  -  REAL  NUMBER 

OUTPUT  PARAMETERS: 

FF  -  CHARACTER  STRING 
NUM  -  INTEGER 
RR  -  REAL  NUMBER 

INPUT  FILES:  NONE 

OUTPUT  FILES:  NONE 

C0^T10N  BLOCKS:  NONE 

ERROR  CONDITIONS:  NONE 

. MAINTENANCE  SECTION . 


DATA  FILES:  NONE 

TEMPORARY  FILES:  NONE 

PRINCIPAL  VARIABLES  AND  ARRAYS: 

SEE  INPUT/OUTPUT  PARAMETERS 
FIRST  -  FLAG  FOR  FIRST  UNIT 
NCHR  -  NUMBER  OF  UNIT  COUNT 
I  DEC  -  LOCATION  IN  NCHR  OF  DECIMAL  TYPE  2 
KNA/KNB  -  ASC! I  NUMBER  RANGE 
KCA/KCB  -  ASCII  LETTER  RANGE 
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C  NSGN  -  (-)  SIGN  FLAG 
C 

C  METHOD: 

C  IDENTIFY  EACH  CHARACTER  BY  TYPE  AND  ACCUMULATE  THEM 
C  ACCORDING  TO  DECODE  TYPE. 

C  PROVIDE  (-)  IF  NOT  CHARACTER  AND  (.)  IF  TYPE  2 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NRFORM 
C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NRFORM 
C 

C  DESCRIPTION; 

C  CONVERT  TO  NORAPS  INPUT  FORM  AND  WRITE  OUT  IN  NORAPS  INPUT  FORMAT 
C 

C  ORIGINAL  PROGRAhWER,  DATE:  HOWARD  LEW  IT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NRFORM 
C 

C  INPUT  PARAMETERS: 

C  ARRAYS  OF  LAT,  LONG 
C 

C  OUTPUT  PARAMETERS : 

C  NONE 
C 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  TAPE  8 
C 

C  COMMON  BLOCKS:  $NRPNT,  $NRAREA 
C 

C  ERROR  CONDITIONS:  BAD  NORAPS  INPUT  WRITE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  SEE  INPUT/OUTPUT  PARAMETERS 
C  MX  =  MAXIMUM  NUMBER  OF  NORAPS  AREAS 

C  ONLY  SET  UP  FOR  AREA  7  12/90 

Q  *mm*********** 

C 

C  METHOD: 

C  CONVERT  LAT,  LONG  TO  NORAPS  INPUT  FORMAT 
C  LAT  NORTH  =  +,  SOUTH  =  - 
0  .ONG  E  =  +,  WEST  =  - 

tia  ~E  NORAPS  INPUT  TO  TAPE  8  FOR  TRANSFER  TO  SAM 
T-aO  OP  OANGES: 

■:::£>>  NRFORM*01  (4  SEPT  91)  —  LEW  IT 
^  A 
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C . END  PROLOGUE . 

C  MAKE  TAPE  13  INPUT  FOR  NORAPS.  INCLUDES  THE  LAT/LONG  OF  THE  NABL 
C  POINTS  AS  WELL  AS  OTHER  NORAPS  CONSTANTS. 

C  CAN  BE  VARIED  BY  NORAPS  AREA  UP  TO  9  AREAS. 

C 

C  EAST  LONGITUDE  =  360  -  LONG 


r-4 1 


NABLS  SDD/IDD  Original 


15  November  1991 


SUBROUTINE  NRLL(AL,NS,AN,EW) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NRLL 
C 

C  DESCRIPTION: 

C  PROVIDE  FOR  INPUT  OF  LAT/LONG  POSITIONS 
C 

C  ORIGINAL  PROGRAhtlER,  DATE:  HOWARD  LEW  IT,  DECEMBER  1990 
C 

C  CONTRACT  NUTBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NRLL(AL,NS,AN,EW) 

C 

C  INPUT  PARAMETERS: 

C  NONE 
C 

C  OUTPUT  PARAMETERS: 

C  AL,NS  =  LAT 
C  AN,EW  =  LONG 
C 

C  INP'JT  F1LES:TAPE1 
C 

C  OUTPUT  FILES:  TAPE2 
C 

C  C0M10N  BLOCKS: 

C 

C  ERROR  CONDITIONS: 

C  NUMEROUS  MESSAGES  TO  USERS  IN  CASE  OF  INVALID  INPUT 


C  UP  TO  NG  ERRORS  ALLOWED  THEN  MESSAGE  TO  GET  HELP  AND  STOP 
C 

C . MAI  NTENANCE  SECT  ION . 

C 


C  DATA  FILES:  NONE 
C 

C  "EMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  LINE-  80  CHARACTER  USER  INPUT 
C  I  TRY  -  ERROR  COUNT 
C 

C  METHOD: 

C  PROMPT  USER  TO  ENTER  FULL  LAT  LONG.  USER  SUBROUTINE  NREXTR  TO 


C  EXTRACT  THE  NUMBERS  AND  CHARACTERS  FROM  THE  ENTRY.  CHECK  THE 
C  ENTRY  FOR  VAL 1 D  VALUES .  VER I FY  W I TH  USER . 

C  IF  NOT  CORRECT  TRY  AGAIN  NG  TIMES. 

C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NROLD 
C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NROLD 
C 

C  DESCRIPTION: 

C  READ  AND  WRITE  POINT  INFORMATION  IN  STANDARD  FORMAT  AS 
C  PERMANENT  FILE  NRNREC(AREA) ,  ID=OP,  SN=SHARDSK 
0 

C  ORIGINAL  PROGRAftlER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NROLD 
C 

C  INPUT  PARAMETERS: 

C  ARRAYS  OF  LAT,  LONG 
C 

C  OUTPUT  PARAMETERS : 

C  NONE 
C 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  PERMANENT  FILE  RECORD  NRNREC(AREA) 

C 

C  COMMON  BLOCKS:  $NRPNT 
C 

C  ERROR  CONDITIONS:  RECORD  NOT  AVAILABLE  CONTINUE  GOING 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  SEE  INPUT/OUTPUT  PARAMETERS 
C 

C  METHOD: 

C  HAVE  THE  USER  INPUT  THE  AREA  TO  BE  UPDATED.  IF  AREA  IS  NOT  DEFINED 
C  FOR  NABL,  NOTIFY  THE  USER  AND  ABORT. 

C  READ  AND  WRITE  ALL  VARIABLES  AS  UNFORMATTED. 

C  NON-CONVERTED  REAL  NUMBERS  OR  AS  INTEGER  REPRESENTATIONS 
C  OF  CHARACTERS.  WRITE  TO  RECORD  NRNREC(AREA)  ON  FILE  ID=OP. 

C  REVERSE  THE  PROCEDURE  WHEN  READING. 

C 

C  RECORD  OF  CHANGES: 

C  <<CHANGE  N0TICE>>  NROLD+C1  (4  SEPT  91)  —  LEWIT 
C  ADD  AREA  8 
C 

C . END  PROLOGUE . 
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SUBROUTINE  NRPROC 
C 

C . START  PROLOGUE . 

C 

0  SUBPROGRAM  NAME:  SUBROUTINE  NRPROC 
C 

C  DESCRIPTION: 

C  PREPARE  PROCEDURE  TO  PASS  AREA  NUMBER  TO  NRNMOD 
C 

C  ORIGINAL  PROGRAM-ER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEC3UENCE): 

C  CALL  NRPROC 
C 

C  INPUT  PARAMETERS: 

C  NONE 
C 

C  OUTPUT  PARAMETERS:  NONE 
C 

0  INPUT  FILES:  NONE 
C 

0  OUTPUT  FILES: 

C  TAPE?  (  PRO) 

0 

C  COrtlON  BLOCKS: 

C  SNRAREA 
C 

C  ERROR  CONDITIONS:  NONE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

0  TEMPORARY  FILES:  NONE 
0 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  NAREA  =  NORAPS  AREA  NUMBER 
0 

C  METHOD: 

0  SET  UP  PROCEDURE  WITH  AREA  NUMBER  INSERTED 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NRRPT 
C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NRRPT 
C 

C  DESCRIPTION: 

C  USER  CHANGE,  DELETE  OR  ADO  POINTS  TO  NRNFIL7 
C 

C  ORIGINAL  PROGRATtlER,  DATE:  HOWARD  LEW IT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NRRPT 
C 

C  INPUT  PARAMETERS:  NONE 
C 

C  OUTPUT  PARAMETERS:  LAT  AND  LONG  IN  ARRAYS  ALAT  AND  ALON 
C 

C  INPUT  FILES:  NONE 
C 

C  OUTPUT  FILES:  TAPE2 
C 

C  COhtlON  BLOCKS : $NRPNT 
C 

C  ERROR  CONDITIONS: 

C  LIMIT  OF  NPT  ENTRY  POINTS,  MESSAGE  AND  CONTINUE 
C  INCORRECT  FORMAT  OF  POINT 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 

U 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  AL  -  USER  INPUT  LAT 
C  AN  -  USER  INPUT  LONG 
C  NS/EW  -  USER  INPUT  N  OR  S  AND  E  OR  W 
C  NMP  -  COUNT  OF  NUMBER  OF  INPUT  POINTS 
C  ALAT  -  LAT  ARRAY 

C  ALON  -  LONG  ARRAY 

C 

C  METHOD: 

C  GET  USER  INPUT  OF  LAT,  LONG  VIA  SUBROUTINE  NRLL 
C  SAVE  MODIFIED  OR  ADDED  DATA 
C  TEST  FOR  ADDITIONAL  INPUT  POINTS  UP  TO  NPT 
C  DELETE  INDICATED  POINT,  CLOSE  SPACE  IN  ARRAYL 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NRSTOP(IP) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  NRSTOP 
C 

C  DESCRIPTION: 

C  TERMINATE  THE  PROGRAM 
C 

C  ORIGINAL  PROGRAWER,  DATE:  HOWARD  LEW  IT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NRSTOP(IP) 

C 

C  INPUT  PARAMETERS: 

C  IP  -  FLAG  1=  NORMAL,  2  =  ABNORMAL 
C 

C  OUTPUT  PARAMETERS:  NONE 
C 

C  INPUT  FILES:  NONE 
C 

C  OLH-PUT  FILES:  NONE 
C 

C  COr-MON  BLOCKS:  NONE 
C 

C  ERROR  CONDITIONS:  NONE 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 

/N 

V/ 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  IP  -  FLAG 
C 

C  METHOD: 

C  STOP  PROGRAM  WITH  MESSAGE  NORMAL  OR  ABNORMAL  TERMINATION 
C  DISCONNECT  TAPE1(  INPUT)  AND  TAPE2(  OUTPUT) 

C  STOP 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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SUBROUTINE  NRYN(*,*) 

C 

C . START  PROLOGUE . 

C 

C  SUBPROGRAM  NAME:  SUBROUTINE  NRYN 
C 

C  DESCRIPTION: 

C  PROVIDE  PROMPT  TO  USER  TO  ENTER  Y  OR  N 
C 

C  ORIGINAL  PROGRATtlER,  DATE:  HOWARD  LEWIT,  DECEMBER  1990 
C 

C  CONTRACT  NUMBER  AND  TITLE: NONE 
C 

C  CLASSIFICATION:  UNCLASS 
C 

C  USAGE  (CALLING  SEQUENCE): 

C  CALL  NRYN(*1,*2) 

C  ALTERNATE  RETURNS: 

C  *1  Y  RESPONSE 

C  *2  N  RESPONSE 

C 

C  INPUT  PARAMETERS:  NONE 
C 

C  OUTPUT  PARAMETERS:  NONE 
C 

C  INPUT  FILES:  TAPE1 
C 

C  OUTPUT  FILES:  TAPE2 
C 

C  COMMON  BLOCKS:  NONE 
C 

C  ERROR  CONDITIONS: 

C  USER  MUST  ANSWER  Y  OR  N,  KEEP  TRYING 
C 

C . MAINTENANCE  SECTION . 

C 

C  DATA  FILES:  NONE 
C 

C  TEMPORARY  FILES:  NONE 
C 

C  PRINCIPAL  VARIABLES  AND  ARRAYS: 

C  YN  -  USER  RESPONSE 
C 

C  METHOD: 

C  PROMPT  USER  TO  ENTER  Y  OR  N.  IF  INCORRECT  KEEP  TRYING 
C 

C  RECORD  OF  CHANGES: 

C 

C . END  PROLOGUE . 
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Appendix  D.  Data  Element  Lists 

Data  Set  N ame  Page 

USER_ONLINE_TERMINAL  D-2 

DATAPOINT_DETAILS/DATAPOINT_NAMELIST  D-3 

NAMELIST_DATA  D-4 

INITIAL_COND_DATA  D-6 

COMPLETE_FORECAST_DATA  D-9 
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Element 

Name _ 

AREA 

ACTION 

P01NT_NR 

LAT_DEG 

LAT_HEMI 

LON_DEG 

LON_HEMI 
LANDORSEA 
Note;  the 
MAX_TAL’ 
Note;  the 


USER  ONLINE  TERMINAL  Data  Elements 


Descr iot ion 

Type/ 

Units 

Lirai ts/Ranee 

Numeric  identifier  of 
larger-scale  model. 

numer i c 

as  displayed 
by  prompt 

Basic  action  -  add, 
delete  or  display  a 
point. 

char  *  1 

"Y"  or  "N" 
response  to 
prompts 

Sequence  number  of 
point  concerned. 

numer i c 

as  displayed 
by  prompt 

Absolute  latitude  of 
data  point . 

real  deg. 

0.0  to  90.0 

Hemisphere  of  lat. 

char ♦ 1 

"N"  or  "S" 

Absolute  longitude 
data  point. 

real  deg. 

0.0  to  180.0 

Hemisphere  of  long. 

cha  r ♦ 1 

"E"  or  "W" 

Does  user  consider  char*5  ’LAND"  or 

point  "wet"  or  "dry".  "OCEAN" 

above  item  is  not  currently  implemented  for  NORAPS. 

Maximum  tau  of  data  to 
be  written  out . 

numer i c 

0  to  L-scal e 
model ’ s  limit 

above  item  is  not  currently  implemented. 
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DATAPOINT_DETAILS/DATAPOINT_NAMELrST  Data  Elements 


El ement 
Name 

Descr iot ion 

Type/ 

Units 

Limits /Ranee 

CAREA 

Identifier  of  larger-scale 
forecast  model . 

char*10 

ex. " NORAPS 7  " 

MAXTAU 

Maximum  hourly  tau  to  be 
output  by  the  model. 

integer 

0  to  L-scal  e 
model ' s  limit 

NRPTS 

Number  of  points  in  file. 

integer 

1-20  (NABL 
dimen ' n  limit) 

PTLAT 

Latitude  of  a  point. 

real 

+/-0.0  -  90.0 

PTLON 

West  longitude  of  point. 

real 

+0.0  -  360.0 

PLAND 

Land  or  Ocean  indicator. 

char *5 

"LAND"  or 
"OCEAN" 

Note 

:  the  three  variables  above 
times  in  each  model-unique 

will  be  repea 
file. 

ted  NRPTS 
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El ement 
Name 
NLTAUF 

NT  ALT 

NFCST 

I  TYPE 

IPT 

JPT 

Note: 

ITAL 

ALAT 

A  LON 
NBLKST 

NR SKEW 

NRTAU 


NAMEL I ST_DATA  Data  Elements 

Type/ 

_ Descr  i  pt  i  on _  L n i  t  s  L  i m i  t  s  /Range 

Final  (highest)  tau  of  NABL  integer  O(dflt)  -  36 

forecast.  (HH)  (dimensioned 

upper  limit) 

Number  output  taus  (number  of  integer  O(dflt)  -  5 
entries  in  ITAL  below).  (dimensioned 

upper  limit) 

Number  of  locations  (number  integer  O(dflt)  -  20 

of  IPT/JPT  or  ALAT/ALON  (dimensioned 

points).  upper  limit) 

Type  of  location  points  pro-  char*2  "lJ"(dflt)  or 

vided  (e.g.  I/J  or  LAT/LON).  "LL” 

Comma  delimited  list  of  1  to  int(20)  larger-scale 

20  I  coordinates  (on  NOGAPS  model  deter- 

or  NORAPS  grid)  where  fore-  mines 

casts  are  desired. 

-  -  -  J  coordinate  equivalent  of  IPT  above  -  -  - 

The  IJ  capability  provided  by  the  three  variables  above 
is  really  DIAGNOSTIC  and  not  required  operationally. 

Comma  delimited  list  of  taus  int(5)  00  -  NLTAUF 

for  the  NABL  forecast  output. 

Comma  delimited  list  of  1  to  real(20)  +/-0.0  to  90.0 
20  latitudes  where  forecasts 
are  desired. 

-  West  longitude  match  for  ALAT  above  -  -0.0  to  360.0 

Comma  delimited  list  of  1  to  integer  five  digits 

20  (WMO)  block-station  idents  (BKSTA) 
corresponding  to  IPT/JPT  or 
ALAT/LON  pairs.  -OPTIONAL- 

Indicates  when  SKEWT  type 
output  from  the  large-scale 
model  input  is  desired.  This 
is  DIAGNOSTIC  INFORMATION 
which  is  written  to  TAPE80. 

Comma  delimited  list  of  1  to  int(8)  00  to  NLTAUF 

8  times  at  which  NRSKEW 
output  (see  above)  desired. 

No  meaning  if  NRSKEW  is  FALSE. 

-  this  table  is  continued  on  next  page  - 


logical  . TRUE .  o  r 

.FALSE. (df 1 t ) 
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E 1 eraent 
N  ame 
NSKEW 


DSKIP 


TSER 


NBSKEW 


FNSKEW 


FNREFC 


Note 


Descr i Dt ion 

Type/ 
Uni  t  s 

Limi ts/Ranee 

Number  of  entries  in  NRTAL 
above.  No  meaning  if  NRSKEW 
is  FALSE. 

integer 

0  (  d  f  1  t  )  -  8 
( d i mens i oned 
upper  limit) 

Indicates  when  limited  data 
is  to  be  written  to  OUTPUT. 
Limits  the  DIAGNOSTIC  infor¬ 
mation  when  set. 

logical 

.TRUE. (df It) 
or  . F AL  SE . 

Indicates  when  DIAGNOSTIC 
time  series  information  is 
to  be  written  to  TAPE70. 
Increases  DIAGNOSTIC  infor¬ 
mation  when  set. 

logical 

.TRUE,  or 
.FALSE. (df It) 

Indicates  when  the  optional 
SKEWT  data  in  NOARL  (metric 
units)  format  is  to  be 
written  to  TAPE? 5. 

logical 

.TRUE,  or 
.FALSE. (df 1 t ) 

Indicates  when  s*andard 

SKEWT  data  in  FNOC  (english 
units)  format  is  to  be 
written  to  TAPE4 2 . 

logical 

.TRUE,  or 
.FALSE. (df 1 t ) 

Indicates  when  standard 
refractive  effects  data  is 
to  be  written  to  TAPE46. 

logical 

. TRUE .  or 
.FALSE. (dfl  t) 

:  PARAMETER  statements  are  correctly  used  to  define 
some  key  array  dimensions  and  loop  control  varia¬ 
bles  used  in  the  NABL  model.  For  examplf’,  the 
upper  dimensions  of  NLTXL'F  and  NFCST  above  aie 
controlled  by  PARAMETER  statements  w  i  t  li  i n  the 
appropriate  program  modules. 
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INITIAL  CON’D  DATA  Data  Elements 


Element  Type/ 

Name  _ Descr  i  pt  ion _  Units  Limits  /Range 


CMOD 

Ident  of  larger-scale  fore¬ 
cast  model . 

char  *6 

"NOGAPS"  or 
"NORAPS" 

CLBL 

Data-set  label.  (UNUSED) 

char  *20 

any  string 

IDIM 

Number  of  geographic  points 
in  this  data  set . 

integer 

1-20  (dimen¬ 
sioned  limit) 

IGMT 

GMT  hour  of  starting  time 
of  the  forecast  (tau=0). 

integer 

0  or  12 

DECL 

Declination  of  sun  at  time 
IFTIM  (below). 

real  deg . 

+/-00.0  to 
+  /-23.5 

lETIM 

GMT  of  tau  zero  (forecast 
ba s  e- t i me ) . 

char *8 

YYMMDDHH 

IFCST 

GMT  of  maximum  forecast  tau. 

char *8 

YYMMDDHH 

NRTAUF 

Max  (last)  tau  in  data  set. 

integer 

0  -  max  f  cs  t 
(NABL  limit  3G) 

PL_IPT 

I  coord,  of  the  data  point 
on  the  forecast  grid. 

real 

grid  limits 

PL_JPT 

J  coord,  of  the  data  point 
on  the  forecast  grid. 

real 

grid  limits 

Note 

:  For  NOGAPS  this  is  the  lower-left  gri 
the  "four  grid  point  set"  containing 
below.  For  NORAPS  this  is  the  nearest 
to  the  LAT/LON  below.  Also,  see  Append 

d  point  in 
the  LAT/LON 
grid  point 
i  X  E . 

PL_LAT 

Latitude  of  the  point. 

real  deg . 

+/-0.0  -  00.0 

PL_LON 

West  longitude  of  the  point. 

real  deg . 

+0.0  -  360.0 

Note 

:  This  is  the  user’s  requested  LAT/LON,  not  the  LAT/LON 
of  I / J  above  (unless  they  coincidently  are  the  same). 

PL_LRS 

"Land-or-Sea"  flag  in  F3 . 1 
(x.y)  format  (see  note  below) 

real 

0.0  to  0.3  or 
1.0  to  1.3 

Note:  "x"  (to  left  of  decimal)  indicates  user's  request  (if 
any)  where  0  is  water  (the  NOGAPS  default)  and  1  is 
land  (the  NORAPS  default).  "y"  (to  right  of  decimal) 
indicates  actual  conditions  at  the  1/J  grid  point  - 
0  is  water,  1  is  land,  2  is  sea-ice  and  3  is  land-ice. 
Also,  see  discussion  in  Appendix  E. 

-  this  table  is  continued  on  the  next  page  - 
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Element 

Name 

Descr i D  t ion 

Type/ 

L'n  i  t  s 

Limi ts/Ranee 

PL_TSFC 

Surface  temperature. 

realdeg.K 

pos i t i ve 

PL_PSFC 

Surface  pressure. 

real (mb. ) 

o  s  i  t  i  V  e 

PL_GWET 

Ground  wetness. 

real  (ND) 

0.0  -  1.0 

PL_TRHT 

Terrain  height. 

real  (m) 

n/a 

PL_DST 

Deep  soil  temperature. 

(Not  avail  from  NOGAPS.) 

realdeg.K 

pos i t i ve 

Note 

:  For  NOGAPS  the  five  variables  above  and  all  of  the 
level  variables  below  are  conditions  interpolated  to 
the  requested  LAT/LON  from  the  "four  grid  point  set" 
containing  that  LAT/LON.  For  NORAPS  they  are  values 
extracted  from  the  single  grid  point  nearest  to  the 
requested  LAT/LON.  Also,  see  discussion  in  Appen¬ 
dix  E . 

TTl 

Level  temperature. 

realdeg.K 

positive 

QQI 

Level  specific  humidity. 

real 
( kg/kg ) 

positive 

uur 

Level  X  component  of  wind. 

real  (mps) 

+/- 

VVI 

Level  Y  component  of  wind. 

real  (mps) 

+  /- 

PPI 

Level  pressure. 

real  (mb . ) 

pos i t i ve 

Note;  the  five  variables  above  are  provided  for  18 
(NOGAPS)  or  21  (NORAPS)  forecast  levels  (from 
lowest  to  highest)  in  the  tau  zero  data  sets  only. 
This  same  tau  zero  information  for  the  first  11 
levels  is  repeated  but  in  a  different  sequence 
following  the  surface  data  below.  This  redundancy 
may  be  worthy  of  elimination  at  some  future  time. 

PSFC  -  same  as  PL_PSFC  above  (one  could  be  eliminated)  - 

TSFC  -  same  as  PL_TSFC  above  (one  could  be  eliminated)  - 

GWET  -  same  as  PL_GWET  above  (one  could  be  eliminated)  - 

ZIL  Surface  roughness  height.  real  (m)  positive 

ID  -  same  as  PL_IPT  above  (one  could  be  eliminated)  - 

JD  -  same  as  PL_JPT  above  (one  could  be  eliminated)  - 


this  table  is  continued  on  the  next  page  - 
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Element  Type/ 

Name  _ Descr  i pt  ion _  Units  L i mi  t s /Range 

TTDN  Exact  tau  of  data  set  which  real  (ND)  positive 

depends  on  time-step  in 
larger-scale  model  which 
might  be  a  few  seconds  off. 


TADB 

-  same 

as 

data 

e 1 ement 

TTl 

above 

QQHR 

-  same 

as 

data 

el ement 

QQI 

above 

SIGP 

-  same 

as 

data 

el ement 

PPI 

above 

UUHR 

-  same 

as 

data 

el ement 

UL'I 

above 

VVHR 

-  same 

as 

data 

e 1 ement 

VVI 

above 

Note:  the  five  variables  above  are  provided  for  the 
first  11  model  (NOGAPS  or  NORAPS)  forecast  levels 
(from  lowest  to  highest)  for  all  tau  data  sets 
including  tau  zero.  See  the  related  note  follow¬ 
ing  data  element  PPl  above. 
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COMPLETE_FORECAST_DATA  Data  Elements 


El ement 
Name 

Descr iot ion 

Type/ 

Units 

L imi t  s/Ranee 

CAREA 

Identifier  of  larger-scale 
forecast  model . 

char  •  1 0 

ex.  "NORAPS7" 

NLTAUF 

Final  (highest)  tau  of  NABL 
forecast . 

integer 

(HH) 

0-36  (dimen¬ 
sioned  limit) 

NTAUF 

Nr.  of  output  taus  (number  of 
entries  in  ITAU  below). 

integer 

O(dflt)  -  5 
(dimen .limit) 

ITAU 

Comma  delimited  list  of  the 
NTAUF  taus  in  NABL  output. 

integer 

00  -  NLTAUF 

NRPTS 

Number  of  forecast  points. 

integer 

1-20  (dimen¬ 
sioned  limit) 

ALAT 

Comma  delimited  list  of  NRPTS 
latitudes  of  the  points. 

real 

+/-0.0  to  90.0 

ALON 

-  West  longitude  match  for  ALAT  above  - 

+0.0  to  360.0 

NRLEVS 

Number  of  levels  in  the 

NABL  forecast. 

integer 

40  (until/ 
unless  chg ' d . ) 

LEVHTS 

Comma  delimited  list  of  NABL 
level  heights  (lowest  first). 

i nt  eger 
(m) 

positive 

TAUSET 

Tau-set  identifier  (based  on 
on  ITAU  sequence.) 

char *5 
(TAUhh) 

ex.  "TAU06", 
"TAU36" 

PTSET 

Point-set  identifier  (based 
on  NRPTS/ALAT/ALON  sequence. 

r  li  a  r  •  5 
( PTnnn ) 

ex.  "PT29  ", 
"PT8",  "PTIOO" 

PL_LRS 

"Land-or-Sea"  flag  in  F3 . 1 
(x.y)  format  (see  note  below) 

r  ea  1 

0.0  to  0.3  or 
1.0  to  1.3 

Note:  "x"  (to  left  of  decimal)  indicates  user's  request  (if 
any)  where  0  is  water  (the  NOGAPS  default)  and  1  is 
land  (the  NORAPS  default).  "y"  (to  right  of  decimal) 
indicates  actual  conditions  at  the  I/J  grid  point  - 
0  is  water,  1  is  land,  2  is  sea-ice  and  3  is  land-ice. 
Also,  see  discussion  in  Appendix  E. 

PL_TRHT  Terrain  height.  real  (m)  n/a 

Note;  the  three  variables  just  above  will  be  repeated 
NRPTS  times  in  each  model-unique  file. 

-  this  table  is  continued  on  the  next  page  - 
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Element  Type/ 


Name 

Descriotion 

U  n  i  t  s 

Limi  t  s/Ranere 

SFCPRES 

Surface  pressure. 

real (mb .  ) 

positive 

GNDWET 

Ground  wetness. 

real 

(ND) 

0.0  -  1.0 

USTAR 

Surface  friction  velocity. 

real 

(mps ) 

pos i t i ve 

THLSTAR 

Surface  scaling  parameter  for 
temperature . 

real 

(K) 

+/- 

MOOBLEN 

Mon in-Obukhov  length  -  scale 
of  boundary  layer  turbulence. 

real 

(m) 

+/- 

SFCRUFF 

Surface  roughness  length. 

real 

(m) 

pos i t i ve 

EVDCTHT 

Evaporative  duct  height. 

real 

(m) 

positive 

SFCMVAL 

Surface  M-refractive  index. 

real 

(ND) 

positive 

Note:  the  eight  variables  above  (starting  with  SFCPRES 
on  the  previous  page  will  be  repeated  NTALT  times 
within  each  point  set  (PTSET). 


OPPRESS 

Grid 

pt.  pressure. 

r  ea  1 

(mb) 

positive 

DDWND 

Gr  id 

pt.  wind  direction. 

int (deg) 

0  -  360 

FFWND 

Gr  i  d 

pt.  wind  velocity. 

r  ea 1 ( mps ) 

positive 

ABSTEMP 

Grid 

pt.  absolute  temp. 

real 

(K) 

positive 

POTTEMP 

Grid 

pt  .  potential  temp. 

r  ea  1 

(K) 

positive 

VPOTEMP 

Grid  pt.  virtual  potential 
t  empr  r  a  t  u  r  e . 

real 

(K) 

positive 

DPT EM P 

Grid 

pt.  dew  point  temp. 

real 

(K) 

pos i t i ve 

SPECHUM 

Grid 

pt.  specific  humidity. 

real (g/kg ) 

positive 

RELHUM 

Grid 

pt.  relative  humidity. 

int. 

(^) 

0  -  100 

MVAL 

Grid 

pt.  M-refractive  index. 

real 

(ND) 

pos i t i ve 

CLOUDS 

Gr  id 

pt.  cloud  fraction. 

real 

(ND) 

0.0  -  1. 

TURBNRG 

Grid  pt.  turbulent  kinetic 
energy. 

real 

(m^/s 

ec^ ) 

pos i t i ve 

LHTFLUX 

Grid 

pt.  latent  heat  flux. 

real 

(wat  ts/m^ ) 

+  /- 

SHTFLUX 

Gr  i  d 

pt.  sensible  heat  flux. 

real 

(wat  t s/m^ ) 

+  /- 
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Note:  the  last  14  variables  on  the  previous  page  (start¬ 
ing  with  OPPRESS)  will  be  repeated  NRLEVS  times  in 
LEVHTS-order  within  each  point-tau  set  (TAUSET). 
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Appendix  E.  NABL  Forecasts  in  the  Coastal  Zone 


The  coastal  zone  presents  particular  problems  to  the  NABL 
forecast  system  as  it  is  now  configured.  The  near-shore  output 
which  NABLS  provides  may  not  be  as  accurate  as  possible  in  some 
cases  and  may  be  unnecessarily  misleading  in  others.  In  this 
appendix  those  problems  are  discussed  and  some  recommended 
changes  to  improve  the  existing  situation  are  presented. 

THE  PROBLEM.  When  a  user  specifies  a  point  ( 1  at i tude/ 1 ong i - 
tude)  near  a  shoreline  at  which  an  NABL  forecast  is  desired,  it 
may  be  assumed  that  the  user  expects  the  base  (the  surface  point) 
of  that  forecast  to  be  specifically  over  water  (the  ocean,  sea, 
gulf,  bay,  estuary  or  whatever)  or  specifically  over  land  (some¬ 
where  above  mean  sea  level).  The  user  might  reasonably  expect 
the  system  to  adjust  the  point  slightly  to  ensure  the  expected 
"wet-"  or  "dry-based"  retrieval;  or,  if  that  is  not  possible,  the 
user  may  reasonably  expect  the  system  to  provide  some  other  type 
of  data  and  to  inform  the  user  of  that  substitution. 

As  the  NABL  System  is  now  implemented,  it  does  not  always 
know  what  the  user  is  expecting;  and,  in  those  cases  when  the 
system  does  know,  the  user  may  receive  something  different  than 
what  was  requested  without  ever  being  notified.  Or,  it  is 
conceivable,  that  the  user  could  be  told  that  a  forecast  is  based 
on  data  over  land  when  in  fact  the  data  better  represents  condi¬ 
tions  over  water  (or  vice  versa).  Such  problems  relate  to; 

a.  the  user  interface  (the  NABL  user  is  not  always 
prompted  to  specify  "land"  or  "sea"  ). 

b.  the  I arger -sea  1  e-mode  I  interface  (the  larger  scale 
models  do  not  require  such  user  input). 

c.  the  1 arger-scal e-mode  1 ’ s  data  retrieval  conventions 
(these  are  inconsistent  and  somewhat  arbitrary  as  will  be  de¬ 
tailed  below). 


d.  the  NABL  output  formats  (NABLS'  failure  to  properly 
display  for  the  user  what  land-sea  information  is  now  being 
provided  by  the  larger  scale  models). 

THE  EXISTING  SYSTEM.  The  NABLS  as  presently  implemented 
prompts  the  user  for  a  "land"  or  "ocean"  assessment  with  respect 
to  the  specified  retrieval  point  only  when  preparing  NAMELIST 
input  for  NOGAPS.  This  is  done  in  program  NGNPNT.  No  such 
prompting  is  done  by  program  NRPPNT  when  preparing  NAMELIST 
input  for  NORAPS.  Nor,  does  either  program  provide  for  a  sophis¬ 
ticated  user  who  may  wish  to  make  a  "sea-ice-"  or  "land-ice- 
based"  retrieval  -  even  though  both  NOGAPS  and  NORAPS  are  capable 
of  making  such  distinctions. 

Neither  NOGAPS  nor  NORAPS  require  a  user's  assessment  as  to 
whether  the  point  is  "wet"  or  "dry"  though  both  are  programmed  to 
accept  such  input.  If  none  is  provided,  NOGAPS  defaults  to 
"ocean"  and  NORAPS  defaults  to  "land"  as  the  "requested  condi- 
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tion".  Both  larg^er  scale  programs  then  proceed  to  ignore  both 
the  request  and  the  default. 

NOGAPS  locates  the  four-grid-point  "square"  containing  the 
requested  latitude  and  longitude  and  interpolates  from  those  four 
corner  grid  points  to  the  requested  geographic  point.  It  returns 
the  interpolated  value  and  the  "base  condition"  of  the  lower-left 
(least  I  and  least  Jl  grid  point  of  the  four-gr  id-point  set. 
This  "retrieved  condition",  as  determined  by  the  model,  may  be 
"water"  (assigned  a  value  of  0),  "land"  (assigned  a  value  of  1), 
"sea  ice"  (assigned  a  value  of  2),  or  "land  ice"  (assigned  a 
value  of  3).  The  "requested  (or  default)  condition"  and  the 
"retrieved  condition"  are  combined  into  a  type  real  " 1 and-or-sea- 
status  flag"  in  x.y  (F3.i)  format  [this  is  variable  PLOC(5)  in 
program  NABL ' s  COMMON  block  /BLANK/]  -  where  x  is  0  or  1  depend¬ 
ing  on  whether  the  user  requested  (or  the  system  defaulted  to) 
"ocean"  or  the  user  requested  "land";  and  y  is  0,  1,  2,  or  3 
(based  on  the  lower-left  corner  as  specified  just  above).  [Note: 
if  NOGAPS  is  provided  I/J  coordinates  rather  than  LAT/LON  coordi¬ 
nate  (as  the  NAMELIST  input  format  permits,  but  the  present  user 
interface  does  not)  there  is  no  interpolation.  The  exact  value 
and  condition  of  the  requested  I/J,  regardless  of  any  user-re- 
quested  condition,  are  extracted  and  returned.] 

Thus,  for  example,  a  user  of  NOGAPS  data  may  request  "land" 
data  and  receive  pure  "ocean"  data  if  the  requested  point  falls 
ever  so  slightly  within  an  all  ocean  grid  point  "square".  In 
such  a  case  the  returned  "land-or-sea-status  (x.y)  flag"  would 
correctly  indicate  the  discrepancy  (would  take  the  value  "1.0"). 
In  another  case,  however,  the  user  might  request  "ocean"  data  and 
receive  data  interpolated  from  three  "land  ice"  points  and  only 
one  "ocean"  point.  In  such  a  case,  and  if  the  "ocean"  grid  point 
were  the  lower-left  point,  the  returned  "land-or-sea-status  flag" 
would  indicate  no  such  problem  (it  would  take  the  value  "0.0"). 

NORAPS  locates  the  grid-point  nearest  to  the  requested 
latitude  and  longitude.  It  returns  the  value  at  that  nearest 
grid  point  regardless  of  whether  or  not  it  is  "wet"  or  "dry".  The 
"retrieved  condition"  of  that  nearest  grid  point  ("land", 
"ocean",  "sea  ice"  or  "land  ice")  is  combined  with  the  "requested 
(or  default)  condition"  and  returned  in  the  same  x.y  (F3.1) 
format  as  used  and  discussed  above  for  NOGAPS. 

Thus,  for  example,  a  user  of  NORAPS  data  may  request  "land" 
data  and  receive  pure  "sea  ice"  data  if  the  requested  point  fall 
ever  so  slightly  closer  to  a  "sea  ice"  point  than  to  a  "land" 
point.  In  such  a  case  the  "land-or-sea-status  (x.y)  flag"  would 
correctly  indicate  the  discrepancy  (would  take  the  value  "1.2"). 

In  none  of  its  output  does  NABLS  take  note  of  and  highlight 
for  the  user  any  discrepancy  between  "requested"  and  "retrieved" 
conditions.  REF_EFFECTS_DATA  and  SKEWT_DATA  outputs  in  FNOC 
format  both  derive  their  "WATER  SURFACE"  or  "LAND  SURFACE"  indi¬ 
cators  from  the  ground  wetness.  This  provides  a  reasonably 
correct  wet  or  dry  data  assessment,  but  fails  to  alert  the  user 
should  the  base  condition  for  the  retrieval  be  different  than 
what  was  requested.  The  SKEWT  output  in  alternate  NOARL  format 
prints  the  real  value  of  the  "land-or-sea-status  flag"  in  F3 . 1 
format  but  does  not  highlight  any  differences  between  the  integer 
and  fractional  components  of  that  flag. 
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RECOMMENDED  CHANGES.  To  improve  NABL  forecasts  in  the 
coastal  zone  the  following  modifications  to  NABLS,  NOGAPS  and 
NORAPS  are  r ecoounended . 

■  Programs  NGNPNT  and  NRPPNT  should  be  modified  (prefera¬ 
bly  combined)  to  provide  a  common  user  interface  which  always 
prompts  for  the  user's  "base  condition"  assessment  as  to  whether 
the  requested  point  location  is  over  "land"  or  "ocean". 

Note;  The  current  version  of  the  NABL  program  is  not  de¬ 
signed  for  and  should  not  be  used  to  prepare  forecasts  over 
sea  ice  or  continental  ice.  However,  two  additional  user- 
input  options  for  "sea-ice"  or  "land-ice"  could  be  useful 
for  a  sophisticated  user  who  recognizes  NABL ' s  limitations 
but  who  may  still  wish  to  take  advantage  of  the  larger-scale 
models'  ability  to  provide  over-ice  data.  Until  such  time 
as  the  Program  NABL  code  has  been  changed  and  certified  for 
higher- 1  at i tude  over-ice  applications,  any  new  user-options 
to  retrieve  sea-  or  land-ice  data  should  be  appropriately 
discussed  in  the  NABL  users  manual  and  the  reader  should  be 
clearly  cautioned  about  their  use. 

•  NOGAPS  should  be  modified  to  always  consider  the  user's 
"base  condition"  assessment  and  to  try  to  retrieve  the  appropri¬ 
ate  type  of  information  without  contamination  from  other  types  of 
"base  conditions".  Thus  the  NOGAPS  program  would  only  interpo¬ 
late  (or  conceivably  extrapolate)  from  points  of  the  user’s 
requested  type.  For  example,  if  two  of  the  four  points  in  the 
four-grid-point  "square"  containing  the  requested  LAT/LON  loca¬ 
tion  were  of  the  requested  type,  data  would  be  retrieved  by 
linear  interpolation  between  those  two  grid  points  only.  A 
similar  scheme  would  be  used  if  only  three  points  matched.  If 
only  one  point  matched,  that  single  point's  values  would  be 
extracted  (or  perhaps  some  extrapolation  using  a  second  nearby 
point  of  the  same  type  could  be  used).  If  all  four  points  in  the 
"square"  were  of  an  unrequested  type,  the  majority  points  only 
(the  two,  three,  or  four  points  of  the  same  unrequested  other 
type)  would  be  used  for  extraction  and  the  type  of  that  majority 
would  be  returned  as  the  fractional  component  in  the 
" 1 and-or-sea-status  flag". 

■  NORAPS  should  be  changed  to  always  extract  values  from 
the  nearest  grid  point  to  the  user's  requested  location  which  is 
of  the  same  type  as  the  user's  stated  "base  condition"  -  un less 
none  of  the  four  grid  points  in  the  f ou r -g r i d-po i n t  "square" 
containing  the  requested  location  are  of  the  requested  type.  In 
the  exception  case  just  stated,  the  absolute  nearest  grid  point 
should  be  used  for  data  retrieval  and  that  point's  type  should  be 
returned  as  the  fractional  component  in  the  " land-or-sea-status 
flag". 

•  For  consistency  with  NOGAPS,  the  NORAPS  "base  condi¬ 
tion"  default  should  be  changed  to  "ocean"  -  which  is  the  more 
prevalent  type  on  a  global  basis  and  more  apt  to  be  of  interest 
to  the  Navy  user  at  sea. 
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■  Program  NABL  should  be  changed  to  boldly  flag  any 
instance  of  a  larger-scale  model  returning  data  based  on  a  grid 
point(s)  of  an  unrequested  type.  This  might  take  the  form: 

*  *  USER  REQUESTED  LAND  DATA  -  MODEL  RETURNED  SEA  ICE  DATA  *  * 

Such  a  statement  would  be  derived  from  the  " land-or-sea~status 
flag"  (variable  PLOC(5)  in  program  NABL's  COMMON  block  /BLANK/) 
and  could  be  routinely  placed  at  the  bottom  of  each  point-tau  set 
of  tabular  data.  If  such  a  statement  were  made  routine,  the 
existing  "LAND  SFC"  or  "WATER  SFC"  output  elements  should  be 
replaced  with  the  actual  ground  wetness  values  ( e . g . "GW=1 . 0000" ) . 
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